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ABSTRACT 


The  safety  provisions  proposed  for  use  in  Canada  for 
limit  states  design  involve  load  factors  to  account  for 
possible  overloads  and  resistance  or  performance  factors  to 
account  for  possible  understrength  of  structural  members. 
The  purpose  of  this  study  was  to  evaluate  the  understrength 
or  <f>  factor  applicable  to  rectangular  tied  reinforced  column 
cross  sections  based  on  a  probabilistic  analysis  of  the 
results  of  a  Monte  Carlo  Study. 

Probability  models  were  described  for  the  major 
variables  affecting  the  cross  sectional  strength.  A  Monte 
Carlo  procedure  was  used  to  develop  a  sample  of  cross 
section  strengths  from  which  the  understrength  factor  was 
calculated.  This  study  showed  that  the  concrete  strength 
variability  and  the  steel  strength  variability  were  the 
major  contributing  factors  to  the  variability  in  cross 
sectional  strength. 

The  understrength  factors  calculated  from  the  results 
of  this  study  were  found  to  be  in  close  agreement  with  the 
understrength  factors  used  in  the  ACI  318-71  Building  Code. 
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CHAPTER  I 


INTRODUCTION 


lil  General 

It  is  generally  recognized  that  there  is  some  degree  of 
uncertainty  in  the  design  equations  used  to  calculate  the 
resistance  of  a  reinforced  concrete  section.  The  strength  of 
a  reinforced  concrete  section  is  calculated  by  the  designer 
as  a  constant  nominal  value  but  it  is  recognized  that  the 
ultimate  strength  of  a  reinforced  column  is  affected  by 
variations  in: 

Concrete  strength 

Steel  strength 

Cross  section  dimensions 

Location  of  steel  reinforcement 

Eccentricity  of  load 

Rate  of  loading 

Amount  of  creep  and  plastic  flow 

The  three  most  common  approaches  that  have  been  used  to 
estimate  the  variability  of  the  ultimate  strength  of  a 
reinforced  concrete  section  are: 

The  technique  of  error  statistics  or  regression 
analysis  applied  to  the  results  of  full  scale  or 
laboratory  tests. 

Direct  statistical  evaluation  of  means  and 
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standard  deviations  from  the  means  and  standard 
deviations  of  the  individual  parameters  involved. 

The  Monte  Carlo  Technique  in  which  the  variables 
affecting  the  cross  section  strength  are  treated  as 
random  variables  and  are  randomly  chosen  and  used  to 
calculate  a  population  of  ultimate  strengths  based  on 
structural  theory. 

The  method  of  error  statistics  has  been  applied  to  sets 
of  test  results  in  various  fields  and  has  been  accepted  as  a 
method  of  analyzing  test  data.  This  method  of  using  test 
data  has  the  disadvantage  of  requiring  many  tests  to  produce 
reliable  results.  More  important,  however,  the  sample  may 
never  be  representative  of  the  population  due  to  testing 
procedures  and  systematic  errors.  Construction  tolerances 
may  not  be  adequately  modeled,  for  example. 

When  the  cross  sectional  strength  can  be  calculated 
with  relatively  simple  analytical  expressions,  standard 
statistical  techniques  can  be  used  to  calculate  the  mean 
error  and  coefficient  of  variation  of  the  cross  sectional 
strength  based  on  the  descriptions  of  the  distributions  of 
the  individual  variables.  This  procedure  becomes  awkward  if 
the  strength  expressions  become  complex. 

The  Monte  Carlo  Technique  has  been  used  to  model  a 
population  of  values  in  various  fields.  This  method  has  the 
disadvantage  of  requiring  a  statistical  description  of  each 
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individual  variable  which  affects  the  final  variable  being 
studied.  The  Monte  Carlo  Technique  has  the  advantage  of 
being  able  to  generate  a  large  size  sample  using  computer 
simulations  rather  than  actual  test  data. 

Since  the  error  statistics  method  of  predicting 
strength  has  been  considered  insufficient  and  too  costly  for 
developing  probability  models  of  cross  section  strength  and 
the  equations  used  to  calculate  the  strength  of  reinforced 
concrete  cross  sections  are  relatively  complex,  the  Monte 
Carlo  Technique  has  become  popular. 

JL.2  The  Monte  Carlo  Technique 

The  Monte  Carlo  Technique  is  a  method  of  obtaining 
information  about  the  total  system  performance  from  the 
individual  component  characteristics.  It  consists  of 
generating  many  total  systems  from  the  component  data  and 
analyzing  the  sample  of  total  systems. 

This  procedure  has  been  used  by  various  researchers  to 
model  the  variability  of  structure  strength  and  loading 
conditions.  Housner  and  Jennings30  have  used  this  procedure 
to  develop  "Artificial  Earthquakes"  from  which  the  various 
effects  of  earthquakes  could  be  measured.  Using  the  data 
generated  with  the  Monte  Carlo  Technique  close  agreement  was 
found  with  actual  measured  values. 

Warner  and  Kabalia72  have  described  a  method  of 
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developing  the  strength  and  serviceability  of  a  real 
structure  using  the  Monte  Carlo  Technigue.  The  strength  of 
an  idealized  axially  loaded  reinforced  concrete  column  was 
calculated  including  the  effects  of  variations  in  the 
material  and  geometric  properties. 

Allen5  has  presented  a  probability  distribution  of  the 
ultimate  moment  and  ductility  ratio  for  reinforced  concrete 
in  bending.  The  ultimate  moment  and  ductility  ratio  were 
obtained  using  prediction  equations  and  probability 
distributions  of  the  parameters.  The  computations  were  based 
on  the  method  of  using  the  Monte  Carlo  Technique  described 
by  Warner  and  Kabalia72.  The  results  showed  that  the 
probability  distributions  of  the  ultimate  moment  and 
ductility  ratio  are  affected  by  material  properties, 
duration  of  loading,  steel  percentage  and  geometric 
properties. 

In  this  study  the  Monte  Carlo  Technigue  was  used  to 
develop  a  probability  model  of  the  strength  of  a  rectangular 
tied  reinforced  concrete  column.  The  actual  probability 
distribution  developed  was  that  for  the  ratio  of  the 
theoretical  load  capacity  to  that  computed  in  accordance 
with  the  ACI  design  equations,  Ptheory/P ACI ,  for  specific 
values  of  e/h  or  eccentricity  of  axial  load.  This  study 
shows  the  effect  of  variations  in  the  concrete  strength, 
steel  strength,  cross  section  dimensions,  location  of 
reinforcing  steel  and  steel  percentage  on  the  probability 
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distribution  of  the  strength  of  a  reinforced  concrete 
section  under  axial  load  and  bending  moment. 

1 . 3  Development  of  the  Under  strength  Factor 

The  ACI  318-71  Building  Code  Requirements  for 
Reinforced  Concrete3  requires  that  the  design  equation 
follow  the  format  of: 

^iYDD+TLL  (1.1) 

Where  <p  is  an  understrength  factor,  R  is  the  nominal 
calculated  resistance  or  strength,  L  and  D  are  the  live  and 
dead  loads  respectively  and  YL  and  YD  are  the  load  factors 
to  account  for  uncertainties  in  the  loads. 

Generally  the  procedure  used  for  determining  the  values 
of  ■  4>  ,  ,  and  yd  has  been  to  rely  on  ’’common  sense  and 
experience”  along  with  a  semi-mathematical  approach.  These 
factors  may  also  be  determined  using  a  logical  mathematical 
approach  using  the  probabilistic  concepts. 

The  first  consistent  proposal  for  design  based  on  the 
concept  of  probability  appears  to  have  been  made  by 
Torroja68.  This  proposal  was  based  on  the  concept  of  limit 
states  in  which  the  design  loads  and  resistance  have  a 
specified  probability  of  being  exceeded. 

Easier10  has  suggested  that  the  coefficient  of 
variation  may  be  used  as  a  probabilistic  but  distribution 
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free  safety  measure.  He  has  proposed  a  rational  scheme  for 
splitting  the  safety  factor  in  the  partial  load  and 
resistance  factors  of  practical  codes.  He  accounts  for  the 
most  uncertain  variables  such  as  workmanship  by  the  use  of  a 
separate  safety  factor  which  may  not  be  defined  explicitly. 
Ang  and  Amin7  have  developed  an  alternative  approach  which 
uses  judgement  factors. 

Cornell*9  has  suggested  a  probability  based  code  format 
known  as  the  first  order  or  second  moment  format.  The  basic 
ACI  Code  design  equation  remains  unchanged  for  this  format, 
whereas  the  code  specified  values  for  loads  and  factors  may 
be  changed  in  absolute  value.  Cornell19  suggests  a  method  of 
calculating  the  values  of  the  loads  and  factors  along  with  a 
method  of  calculating  a  coefficient  of  variation  of  the  in 
place  structure  resistance  and  load  effects. 

The  probabilistic  design  code  suggested  by  Cornell 
implies  that  the  understrength  and  overload  factors  are 
dependent  on  a  predetermined  reliability  factor  and  the 
uncertaintity  in  the  components  which  affect  the  structure 
safety.  In  view  of  this,  equations  must  be  used  to  give  the 
best  estimate  of  the  mean  load  and  strength  conditions  along 
with  their  variability  rather  than  the  current  approach 
based  on  estimates  of  the  material  strengths. 

The  second  moment  theory  developed  by  Cornell  reduces 
to  the  requirement  that  the  mean  safety  margin  be  greater 
than  or  egual  to  a  specified  0  of  its  standard  deviation  so 
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that : 


(1.2) 


Lindas  has  extended  Cornell's  approach  to  code  formats 
of  higher  order  and  demonstrated  a  method  of  calibrating  a 
partial  safety  factor  format  to  Cornell's  as  well  as  Ang  and 
Amin's  format.  Since  it  is  possible  to  choose  6  by 
calibrating  probabilistic  code  formats  to  existing  codes, 
the  parameters  may  be  adjusted  to  yield  designs  comparable 
to  existing  code  designs.  This  leads  to  a  more  acceptable 
implementation  of  probabilistic  code  concepts  initially. 

Siu  et  al.63  have  presented  a  method  of  code 
calibration  which  may  be  used  to  calibrate  probabilistic 
code  formats  with  existing  code  formats  as  well  as  to 
compare  various  probabilistic  code  formats. 

In  this  study  understrength  factors  for  rectangular 
tied  reinforced  concrete  columns  were  calculated  directly 
from  the  distribution  of  column  strength  and  a  probability 
of  understrength  of  1  in  100  and  have  been  compared  with  the 
understrength  factors  calculated  on  the  basis  of  the  first 
order  second  moment  format. 

The  form  of  the  second  moment  format  used  was  that 
developed  by  Cornell,  Lind  and  ACI  Committee  348: 


-BaV, 


R 


(1.3) 


The  derivation  of  this  equation  has  been  reviewed 
Gregor40 . 


by 


CHAPTER  II 


THEORETICAL  BEHAVIOUR  OF  REINFORCED  CONCRETE  SECTIONS 

The  Basic  A ssumptions  for  Analysis 

If  an  analytical  expression  is  to  be  used  to  determine 
the  ultimate  strength  of  a  reinforced  concrete  cross  section 
a  number  of  assumptions  must  first  be  made.  The  following 
basic  assumptions  were  made  for  the  analysis: 

(a)  Plane  sections  remain  plane,  that  is,  the  strain  in 
the  concrete  or  steel  is  directly  proportional  to  the 
distance  from  the  fibre  to  the  neutral  axis. 

(b)  The  concrete  stress  is  a  function  of  the  strain  as 
expressed  by  the  modified  Kent  and  Park  stress  strain  curve 
for  concrete  for  the  theoretical  calculations. 

(c)  The  steel  stress  is  a  function  of  the  strain  as 
expressed  by  an  elastic  plastic  stress  strain  curve. 

W  There  is  no  slip  between  the  concrete  and  steel 
reinforcing. 

(e)  Bending  in  one  plane  is  assumed  and  biaxial  bending 
is  neglected. 

(f)  Stability  failure  of  the  member  is  not  included. 

(g)  The  stiffness  in  bending  of  the  individual  layers 
of  steel  reinforcement  is  neglected. 
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(h)  The  effect  of  duration  of  loading  is  neglected. 

2. 2  The  Stress^ Strain  Relationship  for  Concrete 

The  properties  of  the  compressive  stress  block  of  a 
concrete  flexural  member  may  be  defined  by  the  parameters 
k ,  and  k^  as  shown  in  Figure  2.1.  These  parameters 
depend  on  the  shape  of  the  stress-strain  curve  for  concrete. 

In  North  America  the  most  widely  accepted  stress- strain 

curve  for  concrete  is  that  proposed  by  Hognestad29  which 

consists  of  a  second  order  parabola  up  to  a  maximum  stress 

f”  at  a  strain  e  and  then  a  linear  falling  branch. 

Hognestad 1 s28  curve  was  obtained  from  results  of  tests  on 

eccentrically  loaded  short  columns  in  which  he  found  that 

f"=0.85f • . 
c  c 

There  is  controversy  as  to  whether  the  shape  of  the 
stress-strain  curve  for  concrete  is  affected  by  a  strain 
gradient.  Sturman,  Shah  and  Winter67  concluded  that  the  peak 
occured  at  a  2C%  higher  stress  and  a  50%  higher  strain  for 
eccentrically  loaded  prisms  compared  to  concentrically 
loaded  prisms.  In  Hognestad's  tests  this  was  not  observed. 
There  may  be  no  significant  effect  of  the  presence  of  a 
strain  gradient  but  its  presence,  if  anything,  will  improve 
the  properties  of  the  compression  block.  There  is  no  doubt, 
however,  that  the  presence  of  a  strain  gradient  delays  the 
appearance  of  longitudinal  cracking  in  the  compression  zone. 


' 


1 1 


Figure  2.1 


Compression  Block  Parameters 
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In  columns  the  concrete  is  confined  by  the  ties  to  some 
extent,  depending  on  the  type  of  ties  used.  The  confinement 
due  to  the  ties  does  not  affect  the  concrete  strength  until 
there  has  been  some  yielding  of  the  concrete  to  cause  a  load 
in  the  ties.  At  low  levels  of  stress  the  ties  will  not  be 
stressed  and  therefore  the  concrete  will  act  as  unconfined 
concrete.  Tests  have  shown  that  when  the  stress  in  the 
concrete  approaches  the  maximum  uniaxial  strength, 
deterioration  of  the  concrete  causes  an  outward  expansion 
perpendicular  to  the  load  causing  a  stress  in  the  ties  which 
in  turn  causes  a  confining  pressure.  In  this  case  spiral 
ties  are  more  effective  than  rectangular  ties  since  the 
spiral  is  able  to  exert  pressure  for  its  entire  length 
whereas  the  rectangular  ties  tend  to  exert  pressure  at  the 
corners  and  not  along  their  entire  length.  This  is  due  to 
the  relatively  flexible  bar  between  the  corner  points.  As  a 
result  the  concrete  is  confined  at  the  corners  and  in  the 
centroidal  core  of  the  member.  Even  though  the  rectangular 
ties  are  not  as  effective  as  the  spiral  ties,  they  do 
produce  a  significant  increase  in  ductility  of  the  core  as  a 
whole . 

Some  stress-strain  curves  proposed  for  concrete 
confined  by  rectangular  ties  are  shown  in  Figure  2.2.  In 
Chan* s1 7  trilinear  curve  the  range  OAB  approximates  the 
curve  for  unconfined  concrete  and  the  slope  BC  depends  on 
the  lateral  confinement.  Soliman  and  Yu's64  curve  consists 


Figure  2. 


1  3 


CHAN'S  CURVE 


SOLIMAN  &  YU'S  CURVE 


ROY  &SOZEN'S  CURVE 


Some  Suggested  Stress- Strain  Curves  for 
Confined  Concrete 
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of  a  parabola  and  two  straight  lines.  Values  for  the 
critical  points  are  related  to  the  properties  found  from 
tests  on  eccentrically  loaded  prisms.  Boy  and  Sozen57 
conducted  tests  on  axially  loaded  prisms  and  suggested  that 
the  descending  branch  of  the  stress-st ra in  curve  could  be 
replaced  by  a  straight  line.  The  strain  at  50%  of  the 
maximum  stress  on  the  falling  branch  e^Qc  was  related  to  the 
volumetric  ratio  of  the  transverse  steel. 

Boy  and  Sozen57  concluded  that  rectangular  hoops  did 
not  increase  the  concrete  strength.  Other  investigators  such 
as  Chan17,  Soliman  and  Yu64,  Bertero  and  Felippa11,  and 
Busch  and  Stockl59  have  observed  an  increase  in  strength  due 
to  closely  spaced  rectangular  ties. 

Kent  and  Park36,  on  the  basis  of  experimental  evidence 
have  proposed  the  stress-strain  curve  shown  in  Figure  2.3 
for  confined  and  unconfined  concrete.  This  curve  combines 
many  of  the  features  of  the  previously  described  curves.  The 
ascending  region  AB  is  represented  by  a  second  order 
parabola  in  common  with  the  Hognestad29  curve.  The  confining 
steel  is  assumed  to  have  no  effect  on  the  stress  strain 
relationship  before  the  maximum  stress.  Kent  and  Park35  used 
a  maximum  stress  in  bending  equal  to  f^,,  that  is,  k^=1.0  in 
Figure  2.1.  Sturraan,  Shah  and  Winter* s67  work  suggests  that 
the  value  of  k3=1.C  is  conservative  where  there  is  a  strain 
gradient.  Kent  and  Park36  assume  the  strain,  eQ,  at  maximum 
stress  to  be  0.002  which  is  in  the  range  commonly  accepted 


STRESS 
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STRAIN 


Figure  2.3 


The  Kent  and  Park  Stress-Strai n  Curve  for 
Concrete 
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for  unconfined  concrete.  Confinement  may  increase  the 
maximum  strain  but  this  will  occur  after  the  maximum  stress 
is  reached.  Region  AB  for  the  Kent  and  Park36  curvo  is 
expressed  using: 


f  = 


c 


f' 

c 


(2.1) 


In  this  study  the  value  for  k^  was  taken  as  0.85  based 

on  comparison  with  Hognestad * s29  test  results,  (See  Section 

3.3).  To  allow  compatibility  between  the  ACI  equation  for 

modulus  of  elasticity  and  the  strain  at  about  0.4^,  the 

strain  z  at  a  maximum  stress  of  k  f',  was  taken  as: 
o  3  c 


e 

o 


1.8  f? 
c 


E 

c 


(2.2) 


The  region  of  the  curve  after  the  maximum  stress  is 

linear  from  e  and  f  m  and  is  described  by  the  strain  in 

the  concrete  at  50%  of  the  maximum  stress  as  suggested  by 

Roy  and  Sozen.  The  slope  of  the  falling  branch  increases 

rapidly  with  an  increase  in  concrete  strength.  This  suggests 

that  is  dependent  on  f'.  This  can  easily  be  observed  by 

5ou  c 

the  fact  that  high  strength  concrete  is  more  brittle  than 

low  strength  concrete.  For  concrete  that  is  not  laterally 

restrained,  Kent  and  Park  suggest  that  the  strain  e  at  50% 

5ou 

of  f£  is: 


e  =  3.0  +  0. 002f ' 
5ou  c 

f'  -  1000 
c 


(2.3) 


For  concrete  confined  by  rectangular  ties  the  slope  of 
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the  falling  branch  is  reduced.  This  is  due  mainly  to  the 
restraint  supplied  by  the  ties.  Kent  and  Park35  expressed 
this  in  terms  of  the  ratio  of  the  volume  of  the  ties  to  the 
volume  of  the  concrete  core  within  the  ties.  Kent  and  Park36 
expressed  the  volumetric  ratio  as: 

p"  =  2.0  (b"+dM)  A"  (2.4) 

_ _  s 

b"d"s 


Corley18  suggested  that  the  compression  steel  should  be 
included  in  the  volumetric  ratio.  In  this  study  the 
compression  steel  was  included  in  the  volumetric  ratio  which 
was  expressed  as: 


pM  =  2.0  (b"+d")  A"  +  A'  S 
_ _ s _ s 

b"d"s 


(2.5) 


The  descending  branch  of  Kent  and  Park's36  curve  for 
confined  concrete  may  be  described  by: 


where : 


z 


0.5 


£5oh 


+ 


’5ou 


e 

o 


(2.6) 


(2.7) 


and : 

E50h  =  3/4  *"W  (2’8) 

Kent  and  Park36  assumed  that  confined  concrete  could 
sustain  a  stress  of  0.20^  at  an  infinite  strain  as  shown  by 


the  dashed  line  in  Figure  2.3.  In  this  study  the  descending 
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region  was  assumed  to  continue  to  zero. 

The  tensile  strength  of  concrete  is  usually  neglected 
in  most  flexural  theories  as  well  as  codes  of  practice.  It 
is  reasoned  that  it  may  be  unsafe  to  take  into  account  the 
tensile  strength  of  the  concrete  since  the  concrete  may  be 
cracked  due  to  shrinkage  or  other  reasons  even  before  any 
load  is  applied.  While  the  tensile  strength  of  concrete  is 
small  compared  to  its  compressive  strength  it  has  a  sizeable 
effect  on  the  resistance  and  deformation  of  the  uncracked 
section.  After  the  appearance  of  the  first  cracks  this 
influence  becomes  smaller  and  smaller  as  the  load  increases. 
This  is  due  to  the  fact  that  with  the  advancement  of 
cracking  the  tensile  block  becomes  closer  to  the  neutral 
axis  resulting  in  a  smaller  lever  arm  and  a  negligible 
addition  to  the  moment  capacity. 

In  view  of  the  above  it  was  assumed  that  for  the 
purposes  of  this  study  an  elastic  brittle  stress” strain 
relationship  can  represent  fairly  well  the  behaviour  of 
concrete  in  tension.  An  elastic  brittle  stress- strain 
relationship  can  be  expressed  as  follows: 


for  e  <  £ 

t  -  tr 


(2.9) 


and 


a 


t 


0 


for  e  >  £ 

t  tr 


(2.10) 


The  modulus  of  elasticity  of  concrete  in  tension  was 
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taken  as  the  accepted  value  in  compression: 

E  =  57000  /'T'”' 
c  c 


(2.11) 


The  modulus  of  rupture  was  taken  as  the  accepted  value: 


a  =  7.5  /fr 
tr  c 


(2.12) 


The  complete  stress- strain  curve  for  concrete  used  in 
this  study  is  shown  in  Figure  2.4. 


2.*3  The  Stress-Strain  Relationship  for  Steel 

In  this  study  an  elastic  purely  plastic  stress- strain 
relationship  was  assumed  for  steel  as  shown  in  Figure  2.5. 
The  modulus  of  elasticity  of  steel  was  taken  as  29,000  ksi. 
in  tension  as  well  as  in  compression.  The  steel  stress  was 
assumed  to  increase  to  the  yield  point  and  remain  at  the 
yield  stress  for  any  further  strain.  This  is  a  conservative 
representation  of  the  steel  strength  since  the  effect  of 
strain  hardening  is  neglected. 

2 . 4  Numerical  Method  for  Developing  the  Interaction  Diagram 

The  inter-relationship  between  the  effects  of  the  axial 
load  and  applied  moment  on  a  reinforced  concrete  section  are 
best  shown  by  an  interaction  diagram.  These  diagrams  are  a 
graphical  representation  of  the  envelope  of  the  maximum 
capacities  of  a  reinforced  concrete  section  under  various 
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Figure  2.4 


The  Stress-Strain  Curve  for  Concrete  Used  in 
This  Study 
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Figure  2.5  The  Stress-Strain  Curve  for  Steel  Osed  in  This 

Study 
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axial  load  and  moment  combinations. 

Using  strain  compatibility,  the  moment  curvature 
relationships  were  derived  for  the  section  for  a  number  of 
axial  load  levels  using  the  procedure  described  in  the  next 
few  paragraphs.  The  moment  curvature  relationship  developed 
is  similar  to  that  shown  in  Figure  2.6.  The  maximum  moment 
in  the  moment  curvature  diagram  was  taken  as  the  ultimate 
moment  for  that  given  load.  The  various  values  of  load  and 
ultimate  moment  were  plotted  as  an  interaction  diagram. 

The  calculation  of  the  moment  curvature  diagram  was 

started  by  assuming  a  strain  distribution  across  the  cross 

section  and  determining  the  location  of  the  neutral  axis  and 

the  point  at  which  the  tensile  strains  exceeded  £  .  The 

compression  region  was  then  divided  into  sections  with  equal 

widths  measured  perpendicular  to  the  neutral  axis,  (See 

Figure  2.7).  Using  the  concept  of  linear  strains  in  the 

cross  section  the  strain  at  the  centroid  of  each  section  may 

be  determined.  Fy  assuming  the  strain  is  constant  over  each 

section  the  resulting  stress  and  total  force  over  the  area 

was  determined  with  the  aid  of  the  stress~strain  curve  for 

concrete.  The  total  compressive  force  supplied  by  the 

concrete  may  be  expressed  as: 

ns 

Fc =  E  fcibdx  (2a3> 

i=i 

Assuming  the  maximum  tensile  stress  in  concrete  occurs 

at  a  strain  of  e  and  assuming  a  linear  stress*" strain  curve 

tr 


' 


MOMENT  (k-in. 
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CURVATURE  <£h 


Figure  2.6 


Typical  Moment  curvature  Diagram 


Figure  2.7  Basic  Notation  Used  in  the  Flexural  Analysis  of 

Peinforced  Concrete  Sections 


CROSS  SECTION  STRAINS  STRESSES 
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for  concrete  in  tension,  the  total  tensile  force  may  be 
calculated  using  a  triangular  stress  block.  The  total 
concrete  tensile  force  may  be  expressed  as: 


F  =  c  b  dt 
t  t  — 


(2.14) 


From  the  strain  distribution  the  strain  in  each  steel 

bar  may  be  determined.  Using  the  stre ss- strain  curve  for 

steel  the  stress  in  each  bar  may  be  calculated.  The  total 

steel  force  may  be  expressed  as: 

nb 


F  =  y"  f  .  A. 

st  Z-j  si  1 

i=i 


(2.15) 


The  total  axial  force  resisted  by  the  cross  section  is 
the  algebraic  sum  of  the  concrete  compressive  force,  the 
concrete  tensile  force  and  the  steel  force.  The  total  moment 
that  the  section  is  subjected  to  may  be  determined  by 
summing  the  moments  of  the  above  forces  about  the  centroidal 
axis.  The  moment  may  be  expressed  as: 


m 


ns 

nb 

y  f  .  x.  + 

A-*  ci  l 

i=i 

y  F  x  .  +  F  x  -Pc 

Z-j  st  si  t  t 

i=i 

(2.16) 

where  c  =  the  distance  from  the  tension  steel  to  the 
centroid  of  the  cross  section. 


The  first  three  terms  are  the  moments  of  the  internal 
forces  about  the  tension  steel  and  the  last  term.  Pc,  is  to 
convert  the  moment  to  a  moment  about  the  centroid  of  the 


cross  section 


The  required  points  on  the  axial  load”  moment 
interaction  curve  were  developed  by  selecting  specific  axial 
load  levels  at  which  the  ultimate  moment  was  calculated.  At 
each  axial  load  level  an  initial  strain  at  the  extreme 
compression  fibre  and  initial  curvature  was  assumed.  For  the 
initial  curvature  the  edge  strain  was  incremented  until  the 
sum  of  the  internal  forces  and  the  external  specified  load 
were  balanced  within  a  specified  tolerance.  After  balancing 
the  axial  loads  the  moment  required  to  provide  equilibrium 
was  calculated.  The  curvature  was  then  incremented  and  the 
axial  load  again  balanced  and  the  moment  calculated.  This 
procedure  was  repeated  until  the  maximum  moment  on  the 
moment  curvature  curve  was  calculated. 


The 
gives  a 
accuracy 


subroutine  THEORY  and  flow  diagram  in  Section 
further  description  of  the  above  nrocedure. 
of  this  procedure  is  discussed  in  Section  3.3. 


3.2 

The 


CHAPTER  III 


COMPUTER  PROGRAM  FOR  ANALYSIS 

3 . 1  Description  of  The  Monte  Carlo  Technique 

The  Monte  Carlo  Technique  is  a  method  for  obtaining 
information  about  system  performance  from  the  performance 
data  of  the  individual  components.  This  method  may  be  called 
a  synthetic  or  empirical  method  of  sampling.  It  consists  of 
simulating  many  systems  by  computer  calculation  and  then 
evaluating  the  performance  of  the  overall  system  by 
evaluating  the  performance  of  the  population  of  synthesized 
systems. 

If  a  system  consists  of  many  components  each  with  a 
number  of  values,  a  number  of  systems  could  be  built  to 
measure  the  performance  of  the  system  using  each  component 
value.  Although  this  would  give  an  indication  of  the 
variability  of  the  system,  it  would  generally  be  impractical 
or  uneconomical.  If  there  is  a  relationship  between  the 
total  system  performance  and  each  component  variable,  a 
measurement  of  the  total  system  performance  may  be 
calculated  without  actually  building  the  system.  By  knowing 
the  statistical  properties  of  the  distribution  of  each 
variable  and  drawing  a  value  from  this  distribution  rather 
than  using  measured  values,  it  is  possible  to  calculate  the 
performance  of  a  specified  number  of  synthetic  systems  to 
get  the  variability  of  the  system. 
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This  procedure  is  called  the  Monte  Carlo  Technique  and 
is  shown  graphically  in  the  form  of  a  flow  diagram  in  Figure 
3.1.  The  availability  of  high  speed  computers  has  led  to  the 
popularity  of  this  technique. 

In  this  study  the  Monte  Carlo  Technique  was  used  to 
generate  a  family  of  theoretical  axial  load- moment 
interaction  curves  for  rectangular  column  cross  sections 
using  random  values  of  the  variables  affecting  the  cross 
section  strength.  The  random  value  of  each  variable  was 
based  on  the  statistical  properties  of  each  individual 
variable.  Each  theoretical  curve  was  then  compared  to  the 
ACI  axial  load-moment  interaction  curve  to  obtain  a  sample 
of  ratios  of  the  random  theoretical  capacity  to  that  based 
on  the  ACI  Code,  Ptheory/PACI .  These  ratios  were  eventually 
used  to  calculate  <+>  or  understrength  factors  for  rectangular 
tied  column  cross  sections. 

3^2  Description  of  The  Computer  Program 

The  computer  program  used  in  this  study  is  capable  of 
developing  the  axial  load^moment  interaction  diagram  for 
rectangular  tied  column  cross  sections  with  the  longitudinal 
steel  at  any  location  in  the  cross  section.  The  program  is 
capable  of  developing  the  interaction  diagram  using  the  ACI 
method  and  assumptions  as  well  as  the  theoretical 
interaction  diagram  using  a  theoretical  calculation  of 
strength  based  on  material  and  cross  section  properties. 
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Figure  3 . 1 


The  Monte  Carlo  Technique 
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Figure  3.2  is  a  condensed  flow  diagram  of  the  Monte 
Carlo  program.  The  main  program  consists  of  the  subroutines 
PROP,  ACI,  CURVE,  THMEAN,  RANDOM,  THEORY  and  STAT.  A 
complete  listing  of  the  program  with  its  subroutines  may  be 
found  in  Appendix  D.  Detailed  flow  diagrams  of  the 
subroutines  are  given  in  Appendix  C. 

The  subroutine  PROP  is  used  to  read  and  write  the 
nominal  cross  section  properties.  The  statistical  properties 
of  the  variables  are  read  and  written  in  the  main  program.  A 
complete  description  and  format  of  input  data  is  given  in 
Appendix  E. 

The  subroutine  ACI  is  used  to  calculate  the  ACI  axial 
load^moment  interaction  diagram  using  the  nominal  or 
designer's  values  of  section  and  material  properties.  The 
subroutine  ACI  uses  the  subroutine  ASTEEL  to  calculate  the 
forces  in  the  steel  reinforcement  in  the  cross  section.  The 
capacity  under  pure  axial  load,  balanced  conditions  and  pure 
moment  are  first  calculated.  Using  the  concept  of  linear 
strain  across  the  cross  section  the  axial  load  and 
associated  moment  are  calculated  for  various  strain 
distributions  using  equations  based  on  sections  10.2.1  to 
10.2.5  and  10.2.7  of  ACI  318-713.  Tension  or  compression 
failures  are  classified  by  comparing  the  axial  load  with  the 
axial  load  at  balanced  conditions.  The  value  of  e/h  for  each 
load  level  considered  is  calculated  for  use  in  fitting  a 
curve  to  the  interaction  diagram.  Finally  the  ACI  axial 
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Figure  3.2  Condensed  Flow  Diagram  of  the  Monte  Carlo 

Progra  m 


■ 
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load-moment  interaction  diagram  is  written.  A  condensed  flow 
diagram  of  the  subroutines  ACI  and  ASTEEL  is  shown  in  Figure 

3.3. 

The  subroutine  CDPVE  is  called  a  number  of  times  to  fit 
a  polynomial  curve  to  the  interaction  diagram  developed.  The 
interaction  diagram  is  transformed  into  a  curve  of  axial 
load  vs.  e/h  for  axial  loads  above  the  balance  point  or 
compression  failures  and  a  curve  of  moment  vs.  h/e  for  axial 
loads  below  the  balance  point  or  tension  failures.  The  two 
part  curve  fit  was  used  to  achieve  greater  accuracy  in 
fitting  the  curve  near  the  balance  point.  The  use  of  moment 
rather  than  the  axial  load  was  used  for  the  tension  region 
to  achieve  greater  accuracy  since  e/h  approaches  infinity  as 
P  approaches  zero.  There  was  no  attempt  made  to  force  the 
two  curves  to  coincide  at  the  balance  point  but  the  last 
point  used  for  fitting  the  curve  above  the  balance  point  was 
used  as  the  first  point  for  fitting  the  curve  below  the 
balance  point.  By  using  the  same  point  in  both  curve  fits  a 
close  agreement  was  achieved  at  the  balance  point.  When  the 
subroutine  CDPVE  is  used  to  fit  a  polynomial  to  the 
interaction  diagram  the  calculated  points  with  an  associated 
value  of  e/h  greater  than  3.0  are  eliminated  from  the  curve 
fit  since  these  points  may  cause  large  errors.  Figure  3.4  is 
a  plot  of  the  ACI  interaction  diagram  plotted  from  the  ACI 
calculated  values  and  the  ACI  interaction  diagram  plotted 
from  values  from  the  curve  fit.  Figure  3.4  is  the 
transformed  diagram  with  axial  load  vs.  e/h  and  moment  vs. 
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ACI 


Figure  3.3 


Condensed  Flow  Diagram  of  the  Subroutines  ACI 
and  ASTEEL 


MOMENT  KIP -IN  AXIAL  LOAD  KIPS 
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e/h 


Figure  3.4 


The  ACI  Interaction  Diagram 
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h/e  plotted. 

The  curve  fit  used  for  the  interaction  diagram  required 
a  minimum  of  two  points  above  the  balance  point  and  a 
minimum  of  two  points  between  the  balance  point  and  an  e/h 
value  of  3.0.  This  resulted  in  a  curve  fit  for  three  points 
above  and  below  the  balance  point.  It  was  determined  that  a 

curve  fit  using  a  minimum  of  six  points  resulted  in  a  curve 

\ 

fit  with  virtually  no  error  above  the  balance  point  and  a 
maximum  error  of  about  2.5%  below  the  balance  point  with  the 
general  error  below  the  balance  point  in  the  range  of  1%  or 
less.  On  this  basis  twenty  points  on  the  interaction  diagram 
were  considered  sufficient  to  achieve  a  satisfactory  curve 
fit. 

The  subroutine  CURVE  uses  the  IBM  subroutines  GRATA, 
ORDER,  Ml NV  and  the  modified  IBM  subroutine  MULTR  to  TMULTR. 
These  subroutines  are  described  in  Reference  31.  A  condensed 
flow  diagram  of  the  subroutine  CURVE  is  shown  in  Figure  3.5. 

The  subroutine  THMEAN  uses  the  subroutine  THEORY  to 
calculate  the  theoretical  axial  load-moment  interaction 
diagram  using  the  mean  value  of  the  individual  variables. 
This  subroutine  also  writes  the  interaction  diagram 
calculated.  A  condensed  flow  diagram  of  the  subroutine 
THMEAN  is  shown  in  Figure  3.6. 

The  subroutine  RANDOM  is  a  subroutine  which  combines 
the  IBM  subroutine  GAUSS  and  RANDU  to  calculate  random 


. 


36 


Figure  3.5 


Condensed  Flow  Diagram  of  the  Subroutine  CUFVE 
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Figure  3.6 


Condensed  Flow  Diagram 


of  the  Subroutine  THMEAN 


38 


values  of  each  variable  based  on  the  statistical  properties 
of  each  variable.  These  IBM  subroutines  are  also  described 
in  Reference  31. 

The  subroutine  THEORY  was  developed  to  calculate  the 
theoretical  axial  load~ moment  interaction  diagram  using  the 
subroutines  AXIAL  and  FSTEEL.  A  specific  axial  load  level  is 
chosen  in  THEORY  which  in  turn  calls  AXIAL.  Using  the  axial 
load  level  selected,  a  strain  distribution  is  determined  at 
a  given  curvature  for  which  the  external  load  and  internal 
forces  balance.  For  this  curvature  the  moment  required  to 
develop  the  curvature  is  determined.  The  above  procedure  is 
repeated  with  increasing  curvature  until  the  moment  capacity 
is  determined  at  each  load  level.  This  method  produces  a 
moment  curvature  diagram  similar  to  the  one  shown  in  Figure 
2.6.  The  subroutine  FSTEEL  is  used  by  AXIAL  to  calculate  the 
forces  in  the  reinforcing  steel.  Figures  3.7  through  3.10 
are  condensed  flow  diagrams  of  the  subroutines  THEORY,  AXIAL 
and  FSTEEL.  The  theoretical  interaction  diagram  was  obtained 
as  the  locus  of  the  values  of  Mu  for  each  value  of  P  for 
which  a  moment  curvature  diagram  had  been  computed  as 
explained  in  Section  2.4.  All  comparisons  of  the  theoretical 
strength  with  the  ACI  strength  or  Hognestad*s  tests  were 
done  using  values  of  the  theoretical  strength  after  the 
interaction  diagram  was  subjected  to  a  curve  fit. 

The  subroutine  STAT  is  a  subroutine  used  to  perform  a 
statistical  analysis  on  the  ratio  Ptheory/PACI  for  the 
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Figure  3.7 


Condensed  Flow  Diagram  of  the  Subroutine  THEORY 
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Figure  3.8 


Condensed  Flow  Diagram  of  the  Subroutine  AXIAL 


- 


Condensed  Flow  Diagram  of  the  Subroutine  FSTEEL 


42 


Figure  3.10  Condensed  Flow  Diagram  of  the  Subroutine  STAT 
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various  values  of  e/h  specified.  The  output  from  the 
subroutine  STAT  includes  the  mean,  standard  deviation, 
coefficient  of  variation,  coefficient  of  skewness,  kurtosis, 
minimum  and  maximum  value,  median  and  cumulative  frequency 
table  of  the  ratio  Ptheor y/PACI .  A  condensed  flow  diagram  of 
the  subroutine  STAT  is  shown  in  Figure  3.10. 

3^3  Comparison  of  Theory  With  Test  Results 

The  subroutine  used  to  calculate  the  theoretical  axial 
load-moment  interaction  diagram  was  compared  with  the 
results  of  tests  on  rectangular  tied  columns  reported  by 
Hognestad.  Using  Hognestad’s  column  properties  and  the  total 
eccentricity  of  the  reported  failure  loads  the  mean  ratio  of 
Ptest/Ptheory  was  calculated  to  be  1.0068  with  a  standard 
deviation  of  0.064  when  k^  used  to  define  the  maximum 
compressive  stress,  k  f',  was  taken  equal  to  0.85.  Table  3.1 
is  a  summary  of  the  values  of  the  ratio  Ptest/Ptheory  and 
standard  deviation  for  various  values  of  k^ .  Tables  3.2  and 
3.3  are  the  results  of  a  comparison  of  the  theoretical 
calculations  with  Hognestad * s29  test  results  using  a  value 
of  k3=0.85.  Although  the  lowest  standard  deviation  was 
obtained  for  k3=C.87,  any  increased  accuracy  did  not  warrant 
abandoning  the  traditional  value  of  k^=0.85. 

In  this  study  the  compression  block  was  divided  into 
ten  equal  segments  between  the  extreme  compressive  fibre  and 
the  neutral  axis  with  the  strain  averaged  over  the  segment 
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depth.  A  comparison  of  the  analysis  with  an  analysis  using 
twenty  segments  showed  no  significant  difference  in  the 
ratio  of  Ptest/Ptheory .  The  mean  value  of  P test/Ptheor y  for 
ten  segments  was  1.0068  compared  with  1.005  for  calculations 
using  twenty  segments. 

In  view  of  the  above  the  calculations  in  the  subroutine 
THEORY  were  based  on  ten  segments  with  a  value  of  k^=0.85 
resulting  in  a  mean  value  of  Ptest/Ptheory  of  1.0068  with  a 
standard  deviation  of  0.064.  Any  inaccuracies  due  to  the  use 
of  the  curve  fitting  subroutine  CURVE  are  included  in  these 


statistics. 
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TABLE  3.1 


Comparison  of  Ptest/Ptheory  with  the  value  of 
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Theory  Comparison  With  Hognestad's  Tests  II 
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Theory  Comparison  with  Kognestad's  Tests  III 
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PROBABILITY  MODELS  OF  VARIABLES  AFFECTING  SECTION  STRENGTH 

Concrete  Variability 

iL2.JL.-l  Introduction 

Concrete,  like  all  other  construction  materials,  is 
variable.  This  variability  is  influenced  by  design, 
production  and  testing  procedures.  Research  data  shows  that 
under  current  design  and  construction  techniques  concrete 
which  differs  from  the  specified  strength  is  placed  in 
structures.  These  structures  have  performed  satisfactorily 
due  to  redistribution  of  stresses,  mixing  of  the  under 
strength  concrete  with  over  strength  concrete  within  the 
forms,  and  the  fact  that  the  concrete  strength  increases 
with  age  after  the  time  at  which  tests  are  made.  In  some 
cases  experience  has  lead  to  design  equations  which  result 
in  conservative  designs  even  though  the  assumptions  used  are 
not  entirely  correct. 

The  two  broad  causes  of  variations  in  concrete  strength 
are  variations  in  material  properties  and  variations  in  the 
testing  procedures.  Since  concrete  is  a  heterogeneous 
mixture  of  cement,  water,  coarse  and  fine  aggregate, 
entrained  air,  and  in  some  cases  admixtures,  variations  in 
the  final  concrete  strength  are  inevitable.  Variations  in 
any  one  of  the  ingredients  or  a  combination  of  variations  in 
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more  than  one  ingredient  will  result  in  a  variation  in  the 
final  concrete  strength.  Variation  in  the  water-cement  ratio 
will  cause  significant  strength  variation.  The  water-cement 
ratio  may  be  altered  due  to  poor  control  of  water  content, 
variation  in  moisture  content  or  nonuniformity  of  the 
aggregate.  Variations  in  the  properties  or  proportioning  of 
any  of  the  materials  will  cause  strength  variation.  The 
methods  of  transporting,  placing  and  curing  will  also  affect 
the  final  concrete  strength. 

Variations  in  the  testing  methods  will  lead  to  apparent 
variations  in  the  concrete  strength.  Variations  in  testing 
may  be  due  to  inconsistent  sampling,  nonuniform  fabrication 
of  test  samples  or  poor  handling  and  care  of  fresh  samples 
and  variations  in  temperature  and  moisture  conditions.  Also 
the  preparation  of  the  samples  for  testing  and  the  procedure 
used  in  testing  may  cause  variations  in  the  test  strength. 

The  control  strength  is  affected  by  material  properties 
and  test  procedures  whereas  the  structure  concrete  strength 
will  be  affected  by  the  material  properties  and  placing 
procedures.  This  results  in  different  concrete  strengths  in 
the  test  specimen  and  in  the  structure.  The  concrete 
strength  will  differ  from  place  to  place  in  the  structure 
due  to  different  placing  procedures,  curing  conditions,  and 
the  location  in  the  structure. 
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4.1.2  Distribution  of  Concrete  Strength 

Generally  the  distribution  of  concrete  strength  has 
been  assumed  to  be  a  Gaussian  or  normal  distribution.  ACI 
Committee  2142  found  that  for  practical  concrete  control  the 
normal  distribution  adequately  describes  the  variation  in 
concrete  strength.  Rusch  and  Rackwitz  have  presented  data 
from  an  international  study  of  cube  and  cylinder  tests  which 
also  follows  a  normal  distribution  in  most  cases. 

In  establishing  understrength  factors  for  members  to 
reflect  the  probability  of  the  material  strength  being  lower 
than  the  specified  strength,  the  low  strength  tail  ends  of 
the  curve  are  important.  Eecause  little  data  is  available 
for  these  tail  areas,  the  tail  of  the  curve  must  be 
extrapolated  from  the  central  area  of  the  curve.  The  normal 
distribution  fits  the  data  very  well  for  the  majority  of  the 
data  in  the  central  portion  of  the  curve.  Some  researchers 
have  shown  however,  that  the  normal  curve  does  not  always 
give  the  best  fit  in  the  tail  areas. 

Freudenthal24 ,  Julian33,  and  Shalon  and  Reintz61  have 
shown  that  the  log-normal  distribution  gives  a  better  fit 
for  concrete  strength  in  which  the  control  is  poorer  than 
average  and  should  be  used  where  extreme  values  are 
important.  Shalon  and  Reintz61  have  shown  that  the  normal 
curve  as  a  general  assumption  is  valid  but  in  almost  every 
case  a  skew  towards  the  higher  strengths  was  observed, 
especially  for  cases  of  high  coefficient  of  variation.  Using 
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the  x2  test  as  a  measure  of  discrepancy,  a  discrepancy  was 
observed  between  the  actual  distribution  and  the  normal 
distribution  at  the  5%  level  of  significance  for  a 
coefficient  of  variation  of  23%  whereas  for  a  coefficient  of 
variation  of  14.2%  practically  no  skew  was  observed. 
Freudenthal2 4  suggests  the  use  of  the  log-normal  or  the 
extreme  distribution  to  better  describe  the  tail  area  but 
the  extreme  distribution  has  the  disadvantage  of 
mathematical  complexity. 

Table  4.1  is  a  collection  of  data  from  a  number  of 
statistical  studies  of  concrete  strength.  The  majority  of 
researchers  have  used  a  normal  distribution  due  to  its 
simplicity  and  the  fact  that  in  concrete  control  it  is  the 
central  area  cf  the  curve  that  is  important.  Due  to  this, 
studies  in  concrete  control  are  generally  not  concerned  with 
the  tail  areas  of  the  distribution. 

For  concrete  strengths  with  a  coefficient  of  variation 
of  15%  or  lower  the  normal  curve  describes  the  variation  in 
the  concrete  strength  as  well  as  any  other  distribution.  For 
cases  where  the  coefficient  of  variation  is  greater  than  15% 
a  skewed  distribution  is  observed  for  which  a  log-normal 
transformation  becomes  valid  to  increase  the  accuracy  in  the 
tail  areas  of  the  curve. 
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4.1.3  Statistical  Descript  ion  of  Concrete  Strength 
Variation 

The  average  strength  and  variation  in  strength  of 
concrete  cylinder  tests  may  be  described  by  the  mean, 
standard  deviation  and  coefficient  of  variation.  The 
coefficient  of  variation  has  become  the  accepted  measure  of 
concrete  strength  variation. 

Depending  on  the  control  of  the  concrete  operations  the 
coefficient  of  variation  may  range  from  5%  for  laboratory 
conditions  to  as  high  as  3C%  for  uncontrolled  conditions. 
The  30%  value  is  unacceptable  under  present  construction 
techniques  and  the  5%  value  is  not  practical  for  field 
conditions.  On  the  Skylon  Tower37  at  Niagara  Falls,  Ontario 
coefficients  of  variation  ranging  from  6.8%  to  9.8%  were 
achieved  using  exceptional  control  methods.  This  suggests  a 
minimum  value  for  site  conditions.  The  Bureau  of 
Reclamation2  consistently  achieves  a  coefficient  of 
variation  of  about  15%  which  suggests  a  value  for  better 
than  average  control  or  good  control.  Table  4.1  indicates 
that  the  coefficient  of  variation  in  many  cases  is  between 
15%  and  20%  which  suggests  that  20%  is  a  reasonable  maximum 
value. 

An  ASTK2  task  force  working  on  the  question  of  concrete 
strength  suggested  a  coefficient  of  variation  of  20%  when  no 
control  data  is  available  for  the  average  job.  Figure  4.1 
illustrates  that  the  coefficient  of  variation  varies  but,  on 
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the  average,  is  less  than  20%.  In  this  study  the  levels  of 
control  were  divided  into  three  classes  with  corresponding 
coefficients  of  variation  as  follows: 


Excellent  Control 

10% 

Average  Control 

15% 

Poor  Control 

20% 

The  total  variation  in  concrete  strength  must  include 
the  variation  in  concrete  strength  within  a  single  batch. 
This  in  batch  test  variation  may  be  considered  as  a 
variation  in  testing  procedures  or  a  variation  in  the  actual 
concrete  strength.  The  variation  in  concrete  strength  in  a 
single  batch  will  include  the  effects  of  mixer 
inefficiencies.  Comparison  of  samples  taken  from  different 
locations  in  the  mixer  may  be  used  to  evaluate  the  variation 
within  a  single  batch.  In  this  study  the  levels  of  control 
for  within  batch  tests  were  divided  into  three  classes  with 
corresponding  coefficients  of  variation  as  follows: 


Excellent  Control 

4% 

Average  control 

5% 

Poor  Control 

6% 

Figure  4.1  illustrates  that  the  standard  deviation  and 
the  coefficient  of  variation  are  not  a  constant  for 
different  strength  levels.  Due  to  this  the  mean  strength 
along  with  the  coefficient  of  variation  is  reguired  to 
adequately  describe  the  strength  variation.  The  relationship 
between  the  mean  strength  and  the  standard  deviation  shown 
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in  Figure  4.1  was  developed  using  the  data  from  several 
continents.  The  relationships  shown  by  Murdock,  Erntroy  and 
Rusch47  indicate  that  above  a  certain  value  of  mean  strength 
the  standard  deviation  remains  constant  while  below  this 
value  the  coefficient  of  variation  remains  constant.  The  ACI 
data  indicates  a  constant  standard  deviation  for  all 
strength  levels. 

The  differences  in  values  reported  by  the  different 
sources  may  be  partially  explained  by  the  type  of  data  used. 
The  specimens  of  Erntroy  and  Murdock47  were  6  in.  cubes 
while  the  ACI  specimens  were  standard  cylinders.  The  data 
reported  by  Rusch47  contains  test  specimens  of  both  types. 
Neville46  notes  that  cube  tests  tend  to  be  more  variable 
than  cylinder  tests.  The  relationship  reported  by  Erntroy 
and  Murdock  were  based  on  individual  test  values  whereas  the 
ACI  values  were  based  on  two  specimens  per  test.  The  Rusch 
data  again  contains  both  types  of  data. 

On  the  basis  of  test  data  available  and  reported  it 
appears  that  the  standard  deviation  remains  constant  for 
concrete  strengths  above  a  value  of  3500  to  4000  psi.  and 
the  coefficient  of  variation  is  constant  for  strength  levels 
below  3500  psi. 

4.1.4  Cylinder  Strength  vs.  Design  Strength 


The  average  concrete  strength  required  by  the  ACI 
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Building  Code3  must  exceed  the  value  of  the  design  strength, 
by  at  least: 

400  psi.  if  the  standard  deviation  is  «300  psi. 

550  psi.  if  the  standard  deviation  is  300  to  400  psi. 

700  psi.  if  the  standard  deviation  is  400  to  500  psi. 

900  psi.  if  the  standard  deviation  is  500  to  600  psi. 

If  the  standard  deviation  of  the  test  cylinders  exceeds 
600  psi.  or  if  a  suitable  record  of  test  results  is  not 
available  proportions  shall  be  used  which  provide  an  average 
strength  1200  psi.  greater  than  the  design  strength.  After 
test  data  becomes  available  the  amount  by  which  the  average 
must  exceed  the  design  strength  may  be  reduced  such  that  the 
probability  of  a  test  being  500  psi.  below  the  design 
strength  is  1  in  100  and  the  probability  of  the  average  of 
three  consecutive  tests  being  below  the  design  strength  is  1 
in  100. 

The  amounts  by  which  the  average  strength  must  exceed 
the  design  strength  in  the  ACI  Code  are  based  on  the 
following  criteria: 

1.  The  probability  of  less  than  1  in  10  that  a  random 
individual  strength  test  will  be  below  f^ . 

f  =  f?  +  1.282a 
cr  c 


(4.1) 


. 
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2.  The  probability  of  1  in  100  that  an  average  of  three 
consecutive  strength  tests  will  be  below  f* 


f  =  f'  +  1.343a 
cr  c 


(4.2) 


3.  The  probability  of  1  in  100  that  an  individual 

strength  test  will  be  more  than  500  psi.  below  f*. 

c 

fcr  =  fc  +  2-326a  -500  (4.3) 


where : 

f^  =  the  design  concrete  strength 

fc  =  the  average  cylinder  strength 

a  =  the  standard  deviation  individual  tests 

4.1,5  In-situ  Strength  of  Concrete 

The  concrete  strength  in  a  structure  is  not  clearly 
defined  as  some  specific  multiple  of  the  standard  cured 
cylinder  strength.  Most  researchers  agree  that  the  strength 
of  the  concrete  in  the  structure  is  somewhat  lower  than  the 
standard  test  cylinder  strength. 

Tests  by  Petersons48  on  columns  under  well  controlled 
laboratory  conditions  suggest  that  the  strength  of  the 
concrete  in  the  structure  ranges  from  90%  to  70%  of  the 
standard  cylinder  strength.  Bloem12  suggests  the  strength  of 
the  concrete  in  columns  is  80%  of  the  standard  cylinder 
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strength  for  all  but  the  top  10  in.  of  the  column.  Allen's5 
study  of  beams  failing  in  flexure  suggests  the  strength  of 
concrete  in  the  cases  of  compression  failure  to  be  90%  of 
the  cylinder  strength.  Table  4.2  gives  the  average  ratios  of 
core  strengths  to  cylinder  strengths  from  various 
researchers. 

Petersons46  reviewed  the  data  available  on  core 
strength  as  compared  to  standard  cylinder  strength  and 
concluded  that  the  three  most  important  factors  affecting 
the  strength  of  the  concrete  in  the  structure  are: 

1.  The  strength  level  of  the  concrete”  The  ratio 
between  the  strength  of  the  concrete  in  the  structure  and 
the  standard  cylinder  strength  decreases  as  the  strength 
level  increases. 

2.  The  curing  of  the  concrete”  The  difference  between 
the  minimum  curing  acceptable  and  good  curing  can  be 
approximated  by  a  factor  of  0.9. 

3.  The  location  of  the  concrete  in  the  structure”  Tests 
by  several  researchers  have  indicated  that  the  concrete  in 
the  top  foot  of  columns  is  weaker  than  the  concrete  in  the 
remainder  of  the  column.  This  may  be  explained  by  the 
increased  water  cement  ratio  at  the  top  due  to  water 
migration  after  the  concrete  is  placed.  The  reduction  in 
strength  is  of  the  order  of  15%  of  the  strength  of  the 
remainder  of  the  column. 


Table  4.2 


Concrete  Strength  in  Structure  vs.  Cylinder  Strength 


Researcher 


Core  Strength 

Ratio  - - 

Cylinder  Strength 


|  Kaplan 

— i - 

1 

1 

0.74 

8 

S 

1 

0.96 

1 

\ 

a 

0.90 

1  Petersons 

1 

1 

i 

0.90 

1 

1 

| 

0.88 

|  Bloem 

1 

1 

a 

0.83 

|  Campbell  and  Tobin 

1 

! 

1 

0.87 

i 

I 

(1 

1 

61 


The  reduction  in  the  concrete  strength  in  the  structure 
is  partially  offset  by  the  requirement  that  the  average 
cylinder  strength  must  be  from  700  to  900  psi.  greater  than 
the  design  strength  to  meet  existing  design  codes.  Based  on 
this  observation  and  on  the  equations  from  Allen  and  Bloem40 
the  mean  strength  for  minimum  acceptable  curing  may  be 
expressed  as: 

f  ,  „  ,  s  =  (0. 675f '  +  1.1)  ksi  (a  a) 
c(structure)  c 


h . 1 . 6  Probability  Model  for  Concrete  Strength 


The  variation  in  concrete  strength  was  described  with  a 


normal  distribution  and  a  mean  value  of: 


with  a  coefficient  of  variation: 


V 


2  2 
V2  +  V3 


(4.5) 


(4.6) 


where : 

V  =  the  variation  in  test  cylinder  strengths 

1 

V  =  the  variation  between  real  strength  of  cylinders 

2 

and  measured  cylinder  strengths,  "in°test 
variation" 


V  =  the  variation  of  in*»situ  strength  relative  to 
3  ' 

cylinder  strength 


' 
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The  basic  cylinder  strength  variation  was  taken  as  0.15 
with  a  basic  in-test  variation  of  0.04  and  a  variation  of 
0.10  for  differences  between  in-situ  and  test  cylinder 
strengths.  Checks  were  also  made  with  test  cylinder 
variation  of  0.10  and  0.20. 

4 . 2  Reinforcing  Steel  Variability 

The  variability  of  the  strength  of  the  reinforcing 
steel  was  described  with  a  normal  distribution  as  well  as  a 
modified  log-normal  distribution.  The  complete  discription 
of  the  reinforcing  steel  strength  distribution  used  is  given 
in  Appendix  A. 

4^3  Cross  Section  Dimensional  Variability 

Introduction 

Geometric  imperfections  are  the  variations  in  the 
dimensions,  shape,  lines,  grades  and  surfaces  of  as-built 
structures  compared  to  the  specified  values.  Variations  in 
cross  section  dimensions,  verticality  of  columns, 
misalignment  and  intial  curvature  of  columns  are  inevitable 
in  structures.  Geometric  imperfections  arrise  during  each 
phase  of  the  construction  process.  Variations  in  the  size 
and  shape  are  particularly  dependent  on  the  size,  shape  and 
quality  of  the  forms  used  for  manufacture.  Setting  out  and 
assembly  affect  the  geometry  of  the  overall  structure  and 
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are  dependent  on  construction  techniques  and  construction 
and  inspection  personnel. 

Data  from  field  measurements  of  imperfections  is  needed 
for  various  purposes  such  as  for  the  evaluation  of  specified 
tolerances,  construction  performance  and  theoretical 
probability  models.  It  is  important  that  data  be  collected 
which  is  complete  and  uniform.  Unfortunately  at  present 
there  is  not  a  uniform  method  of  collecting  and  reporting 
this  data.  Without  some  degree  of  standardization  it  is 
difficult  to  compare  the  results  of  measurements  made  by 
various  investigators  with  any  degree  of  reliability. 

4.3.2  Probability  Model  for  Cross  Section  Dimensions 

The  variation  in  column  cross  section  dimensions  has 
been  reported  by  Tso  and  Zelman69.  Their  results  are 
summarized  in  the  histogram  in  Figure  4.2.  The  dimensional 
measurements  were  made  to  the  nearest  1/4  in.  in  conjunction 
with  a  study  of  the  strength  variation  in  concrete.  The  data 

is  based  on  299  columns  from  8  buildings.  The  nominal 

« 

dimensions  ranged  from  12  in.  to  30  in..  Usually  two 
measurements  were  made  at  each  of  five  levels  over  the 
storey  height  of  the  column.  The  mean  variation  was  found  to 
be  +  0.06  in.,  that  is,  the  width  or  thickness  averaged  0.06 
in.  larger  than  the  specified  value  with  a  standard 
deviation  of  C.28  in..  Tso  and  Zelman*s69  measurements 
indicate  the  distribution  of  dimensional  variations  to  be 
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Histogram  of  Cross  Section  Dimensional 
Variation  Reported  by  Tso  and  Zelman 


Figure  4.2 
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normal . 

The  variation  in  the  dimensions  of  one  size  of  column 
has  been  reported  by  Hernandez  and  Martinez28.  Their  results 
are  summarized  in  the  histogram  in  Figure  4.3.  The 
measurements  were  made  at  five  levels  over  the  storey  height 
of  the  column.  At  each  level  the  width  and  thickness  were 
measured  at  each  face  and  at  the  centre  line  of  the  column, 
seventeen  columns  were  studied  with  a  nominal  cross  section 
of  11.811  in.  (30  cm.)  by  19.685  in.  (50  cm.).  A  mean 
variation  was  found  to  be  +  0.15  in.,  that  is,  the  width  or 
thickness  averaged  0.15  in.  larger  than  the  specified  value 
with  a  standard  deviation  of  0.157  in..  A  normal 
distribution  also  describes  the  variation  in  column 
dimensions  reported  by  Hernandez  and  Martinez28. 

Fiorato23  has  reported  a  mean  deviation  of  0.0118  in. 
(0.3  mm.)  to  0.276  in.  (7.0  mm.)  with  a  standard  deviation 
ranging  from  0.063  in.  (1.6  mm.)  to  0.154  in.  (3.9  mm.)  for 
precast  beams  and  columns  ranging  in  size  from  7.87  in.  (200 
ram.)  to  23.62  in.  (600  mm.).  These  values  are  based  on  a 
collection  and  comparison  of  published  data  from  field 
measurements,  primarily  from  Sweden.  These  values  may  not  be 
considered  comprehensive  but  do  give  an  indication  of  the 
variations  which  may  occur  in  prefabricated  structures. 

As  stated  earlier,  due  to  inconsistencies  in  measuring 
and  reporting  techniques,  comparison  of  data  on  column 
dimensions  from  various  researchers  is  difficult.  The 
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Histograms  of  Cross  Section  Dimensional 
Variation  Reported  by  Hernandez  and  Martinez 
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majority  of  researchers  are  interested  in  actual 
construction  tolerances  in  which  a  maximum  and  minimum 
tolerance  is  reported  rather  than  a  mean  value  and  its 
standard  deviation. 

In  this  study  a  normal  distribution  was  used  to 
describe  the  variation  in  column  dimensions  with  a  mean 
value  of  +  0.06  in.  and  a  standard  deviation  of  0.28  in.. 
Tso  and  Zelman's69  results  were  used  since  they  are  based  on 
North  American  data  and  are  based  on  a  larger  sample  size 
than  that  obtained  by  Hernandez  and  Martinez28. 

4 . 4  Reinforcing  Steel  Placement  Variability 

Redkop53  has  developed  models  describing  the  error  in 
placing  reinforcing  steel  in  rectangular  tied  columns  based 
on  test  data  from  measurements  on  several  columns  in  several 
buildings.  He  describes  the  variation  in  steel  placement 
with  respect  to  the  specified  cover  for  the  steel  in  the 
exterior  layers  and  the  specified  position  for  the  interior 
steel.  The  error  in  steel  placement  may  be  described  by  the 
normal  probability  distribution.  Redkop53  observed  that  the 
placement  error  was  a  function  of  the  column  size  as  well  as 
construction  practices.  Since  statistical  data  available 
does  not  suggest  a  complicated  relationship,  a  linear 
relationship  between  column  size  and  placement  error  was 
assumed  with  a  normal  distribution  of  scatter. 


Based  on  Redkop1 s53  data  the  error  in  placement  of  the 
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interior  steel  may  be  described  with: 

e  =  0.04  in. 
n 

a  =  0.2035  +  0.0329  h  (4.7) 

The  placement  of  the  steel  in  the  exterior  layers  may 
be  described  with: 

C  =  C  +  0.250  +  0.0039  h  (4.8) 

a  sp 

a  =  0.166 

where : 

eR  =  placement  error  of  interior  steel  in  inches. 

°  =  standard  deviation  in  inches. 

h  =  column  dimension  perpendicular  to  the  neutral 
axis. 

c_  =  actual  cover  of  exterior  steel  in  inches. 

CL 

c  =  specified  cover  of  exterior  steel  in  inches. 

S  p 

Based  on  Bedkop's53  data  the  mean  variation  in  concrete 
cover  of  the  exterior  steel  is  *  0.315  in.,  that  is,  the 
actual  cover  on  the  average  is  0.315  in.  larger  than  the 
specified  cover,  with  a  standard  deviation  of  0.166  in.. 
Hernandez  and  Martinez  report  a  mean  variaton  of  +  0.473  in. 
with  a  standard  deviation  of  0.13  in..  The  smaller  standard 
deviation  of  the  Mexican  data  is  due  to  measurements  being 
taken  only  from  one  size  of  column  whereas  the  measurements 
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reported  by  RedkopS3  were  taken  from  various  sizes  of 
columns.  Figure  4.4  is  a  histogram  summarizing  the  results 
of  concrete  cover  reported  by  Hernandez  and  Martinez17.  The 
normal  distribution  may  be  used  to  describe  the  variation  in 
the  concrete  cover  for  both  sets  of  measurements. 

In  this  study  the  error  in  steel  placement  was 
described  by  Equations  4.7  and  4.8  with  a  normal 
distribution.  Negative  cover  is  not  a  problem  since  with  1 
1/2  in.  nominal  cover  negative  cover  does  not  occur  before 
the  value  of  cover  is  the  mean  value  minus  10.54  standard 
deviations.  The  probability  of  the  value  of  cover  being  the 
mean  value  minus  eight  standard  deviations  is  approximately 
6.22x10*6. 
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COVER  (cm) 


Figure  4.4 


Histogram  of  Variation  in  Concrete  Cover 
Reported  by  Hernandez  and  Martinez 


CHAPTER  V 


THE  MONTE  CARLO  STUDY 

5 . 1  Size  of  Columns  and  Reinforcement  Studied 

For  this  study  the  size  of  columns  and  reinforcement 
selected  was  based  on  a  limited  study  of  columns  in  existing 
structures.  A  column  take-off  was  done  on  five  buildings 
including  a  high  rise  office  building,  a  parking  garage,  a 
university  building,  a  hospital  and  an  industrial  type 
building. 

Figure  5.1  is  a  histogram  of  the  frequency  of  column 
size  vs.  column  size.  This  histogram  indicates  that  the 
majority  of  columns  are  24  in.  or  smaller.  The  high 
percentage  of  columns  in  the  52  in.  to  56  in.  range  is  due 
to  the  small  number  of  buildings  studied  in  which  one  was  a 
high  rise  with  large  columns  throughout.  From  the  histogram 
of  column  sizes  the  12  in.  column  and  the  24  in.  column  were 
taken  as  representative  of  the  smaller  and  larger  sizes  of 
columns. 

Representative  reinforcing  steel  percentages  were 
chosen  in  the  same  manner  as  the  column  sizes.  Figure  5.2  is 
a  histogram  of  the  reinforcing  steel  percentage  used  in  all 
columns.  Figures  5.3  through  5.5  are  histograms  of  steel 
percentages  used  in  the  various  sizes  of  columns.  From  these 
histograms  it  can  be  seen  that  the  most  commonly  used  steel 
percentage  ranges  from  1%  to  3%»  Based  on  these  histograms  a 
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Figure  5 . 1 


Histogram  of  the  Frequency  of  Column  Sizes 
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%  STEEL 


Figure  5.2 


Histogram  of  the  Percentage  of  Feinforcing 
Steel  in  All  Columns 


7  4 


Figure  5.3 


Pistoqram  of  the  Percentage  of  Reinforcing 

Steel  in  Columns  Less  Than  16  in. 


75 


0 


1.5 


3.0  4.5 

%  STEEL 


6.0 


7.5 


Figure  5.4  Histogram  of  the  Percentage  of  Feinforcing 

Steel  in  Columns  16  in.  to  24  in. 
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Figure  5.5 


Histogram  of  the 
Steel  in  Columns 


Percentage  of  Reinforcing 
24  in.  to  36  in. 
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steel  percentage  of  1%  was  chosen  for  a  lower  limit  and  a 
steel  percentage  of  3%  was  chosen  for  an  upper  limit. 

The  final  column  cross  sections  chosen  are  shown  in 
Figure  5.6.  The  basic  column  was  12  in.  by  12  in.  with  a 
nominal  steel  percentage  of  1%.  A  24  in.  by  24  in.  column 
was  chosen  to  have  a  low  variability  of  strength  with  a 
nominal  steel  percentage  of  3%.  The  nominal  or  designer's 
concrete  and  steel  strengths  were  3000  psi.  and  40000  psi. 
respectively.  These  strengths  and  properties  were  chosen  to 
get  an  estimate  of  the  upper  and  lower  bounds  of  the 
variabilities. 

Interaction  diagrams  for  the  two  sections  are  presented 
in  Figures  5.11  and  5.12  and  will  be  dicussed  more  fully  in 
Section  5.4.  The  balanced  eccentricity,  e^/h,  was  0.4  for 
the  12  in.  column  and  0.5  for  the  24  in.  column.  The  columns 
are  fully  described  in  Appendix  B  with  their  nominal 
properties  and  the  mean  values  and  standard  deviations  of 
the  variables  affecting  column  strength. 

5^2  Size  of  Sample  Studied 

For  this  study  a  sample  size  was  required  which  would 
give  reasonable  results  compared  to  a  large  sample  size 
without  using  an  excessive  amount  of  computer  time.  Sample 
sizes  of  1000,  50C  and  200  were  used  to  determine  the 

smallest  practical  sample  size. 
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Figure  5.6 


Final  Column  Cross  Sections  Studied 
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Figure  5.7  is  a  plot  of  the  mean  value  of  Ptheory/PACI 
vs.  e/h  for  each  sample  size.  The  mean  value  is  practically 
independent  of  the  sample  size  so  that  any  of  the  sample 
sizes  could  be  used  to  determine  the  mean  value  of 
Ptheory/PACI. 

Figure  5.8  is  a  plot  of  coefficient  of  variation  of  the 
ratio  Ptheory/PACI  vs.  e/h  for  the  three  sample  sizes.  The 
coefficient  of  variation  for  the  sample  size  of  500  is 
practically  the  same  as  the  coefficient  of  variation  for  a 
sample  size  of  1C00  over  the  range  of  e/h  less  than  1.0. 
Since  a  good  correlation  was  found  between  the  sample  size 
of  500  and  the  sample  size  of  1000  below  an  e/h  value  of  1.0 
the  sample  size  of  500  was  acceptable  when  the  mean  and 
coefficient  of  variation  were  needed  as  output. 

Figure  5.9  is  a  plot  of  the  coefficient  of  skewness  of 
the  ratio  of  Ptheory/PACI  vs.  e/h  for  the  three  sample 
sizes.  The  coefficient  of  skewness  for  a  sample  size  of  500 
is  not  significantly  different  from  the  coefficient  of 
skewness  for  a  sample  size  of  1000. 

Tables  5.1  through  5.4  are  tables  of  comparison  of  the 
mean  values,  coefficients  of  variation,  coefficients  of 
skewness  and  kurtosis  of  the  ratio  of  Ptheory/PACI  for 
sample  sizes  of  200,  500  and  1000. 

All  the  tables  and  graphs  of  comparison  indicate  no 
meaningful  increase  in  accuracy  in  using  a  sample  size  of 
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Mean  Value  of  the  Ratio  Ptheory/PACI  vs.  e/h 
for  Sample  Sizes  of  200,  500  and  1000  for  a  12 
in.  Sguare  Column  and  Modified  Log-normal  Steel 
Strength  Distribution 
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Figure  5.8  Coefficient  of  Variation  of  the  Patio 

Ptheory/PACI  vs.  e/h  for  Sample  Sizes  of  200, 
500  and  1000  for  a  12  in.  Square  column  and 
Modified  Log-normal  Steel  Strength  Distribution 


COEFFICIENT  OF  SKEWNESS 


82 
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Coefficient  of  Skewness  of  the  Patio 
Ptheory/PACI  vs.  e/h  for  Sample  Sizes  of  2C0, 
500  and  1000  for  a  12  in.  Square  Column  and 
Modified  Log-normal  Steel  Strength  Distribution 


Figure  5.9 
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Table  5 . 1 


Comparison  of  the  Mean  Value  of  the  Ratio  Ptheory/PACI 
Sample  Sizes  of  290,  500  and  1000 
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Table  5.2 


Comparison  of  the  Coefficient  of  Variation  of  the  Patio 
Ptheory/PACI  for  Sample  Sizes  of  200r  500  and  1000 
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Table  5 . 3 


Comparison  of  the  Coefficient  of  Skewness  of  the  Patio 
Ptheory/PACI  for  Sample  Sizes  of  200,  500  and  1000 
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Normal  and  log-normal  distributions  have  coefficients  of 
skewness  of  0.0  and  0.5  to  1.5  respectively. 
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Table  5.4 


Comparison  of 
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of 
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A  normal  distribution  has  a  Kurtosis  of  3.0 
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1000  over  a  sample  size  of  500.  On  this  basis  sample  size  of 
500  was  used  for  all  subsequent  calculations. 

5.3  Results  of  The  Monte  Carlo  Simulation 


5.3.1  General 

In  this  Monte  Carlo  study  the  relationship  between  the 
theoretical  axial  load-moment  interaction  diagram  and  the 
ACI  axial  load-moment  interaction  diagram  was  determined. 
This  relationship  was  calculated  using  the  Monte  Carlo 
Technique  to  give  the  mean  ratio  of  the  theoretical  strength 
divided  by  the  PCI  strength  along  with  its  standard 
deviation  at  various  e/h  values.  From  the  mean  ratio,  the 
standard  deviation  and  the  type  of  distribution  an 
understrength  or  <{>  factor  was  calculated. 

To  aid  in  the  development  of  an  unde rstrength  factor 
the  effect  of  the  variation  in  concrete  strength,  steel 
strength,  cross  section  dimensions  and  location  of  the 
reinforcing  steel  was  studied.  In  addition  the  effect  of  the 
type  of  distribution  of  steel  strength  used  was  studied 
using  a  normal  distribution  and  a  modified  log-normal 
distribution . 


5.3.2.  The  Effect  of  Steel  Strength  Distribution  Used 


The  effect  on  the  strength  of  the  column  cross  section 
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of  the  type  of  distribution  of  steel  strength  was  studied 
using  a  normal  and  a  log-normal  distribution  of  steel 
strength.  Both  types  of  distribution  can  be  fitted  to  the 
data  on  steel  strength  as  shown  in  Appendix  A.  Tables  5.5 
through  5.8  are  tables  of  comparison  of  the  mean  values, 
coefficients  of  variation,  coefficients  of  skewness  and 
measure  of  kurtosis  of  the  ratio  Ptheory/PACI  for 
calculations  based  on  steel  strength  normally  distributed 
and  steel  strength  which  follows  a  modified  log-normal 
distribution. 

The  distribution  assumed  for  the  variation  in  the  steel 
strength  did  not  significantly  affect  the  mean  ratio  of  the 
theoretical  strength  to  the  ACI  strength  as  shown  in  Table 
5.5  but  did  affect  the  distribution  of  the  ratio  in  the 
tension  failure  region  of  the  interaction  diagram.  When  the 
modified  log-normal  steel  strength  distribution  was  used, 
the  distribution  of  the  ratio  of  theoretical  strength  to  the 
ACI  strength  approached  a  log-normal  distribution  for  values 
of  axial  loads  below  the  balance  point.  If  a  normal 
distribution  of  the  steel  strength  were  used,  this  ratio  was 
normally  distributed.  This  is  shown  by  the  coefficient  of 
skewness  given  in  Table  5.7.  For  normally  distributed  steel 
yield  strengths  the  coefficient  of  skewness  remained  close 
to  zero  throughout  corresponding  to  a  normal  distribution. 
With  the  log-normal  assumption  the  coefficient  of  skewness 
approaches  1.0  for  tension  failures  corresponding  to  a  log¬ 
normal  distribution,  (  See  also  Figure  5.9  ).  The  use  of  the 
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Table  5.5 


Comparison  of  the  Mean  Value  of  the  Patio  Ptheory/PACI  for  a 
Normal  and  a  Modified  Log-normal  Steel  Strength  Distribution 
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Table  5.6 


Comparison  of  the  Coefficient  of  Variation  of  the  Patio 
Ptheory/PACI  for  a  Normal  and  a  Modified  Loq-normal  Steel 
Strength  Distribution 
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Table  5.7 


Comparison  of  the  Coefficient  of  Skewness  of  the  Fatio 
Ptheory/PACI  for  a  Normal  and  a  Modified  Log-normal  Steel 
Strength  Distribution 
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Table  5.8 

Comparison  of  the  Measure  of  Kurtosis  of  the  Patio 
Ptheory/PACI  for  a  Normal  and  a  Modified  Log-normal 
Distribution  of  Steel  Strength 
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modified  log-normal  steel  strength  distribution  resulted  in 
a  larger  4>  factor  at  the  1%  level  of  probability  of  failure 
than  that  for  the  normal  distribution  of  steel  strength  in 
the  tension  region  of  the  interaction  diagram.  The 
calculation  of  this  term  is  discussed  in  Section  5.5. 

The  type  of  steel  strength  distribution  used  did  not 
significantly  affect  the  calculated  value  of  the  <j>  factor  at 
the  1%  level  of  probability  of  failure  in  the  compression 
failure  region  of  the  interaction  diagram.  The  distribution 
of  the  ratio  of  the  theoretical  strength  to  the  ACI  strength 
in  the  compression  failure  region  followed  a  normal 
distribution  for  both  types  of  steel  strength  distribution. 
This  may  be  explained  by  the  failure  in  the  compression 
region  being  dependent  on  the  concrete  strength  rather  than 
the  steel  strength. 


5.3.3  The  Effect  of  The  Concrete  Strength  Variation 

The  effect  of  the  coefficient  of  variation  of  the 
concrete  cylinder  strength  was  studied  by  keeping  all  other 
variables  at  their  mean  values  and  using  a  cylinder  strength 
variation  of  10%,  15%  and  20%.  The  overall  coefficients  of 
variation  of  in=situ  strength  were  13.6%,  17.6%  and  22%  as 
computed  per  Egn.  4.6  in  Section  4.1.o.  Tables  5.9  through 
5.11  are  tables  of  comparison  the  mean  values,  coefficients 
of  variation  and  skewness  of  the  ratio  Ptheory/PACI  for 
various  values  of  e/h  and  cylinder  strength  coefficients  of 
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Table  5.9 


Comparison  of  the  Mean  Value  of  the  Patio  Ptheory/PACI  for 
Concrete  Cylinder  Strength  Coefficients  of  Variation  of  10%, 
15%  and  20% 
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Table  5.10 


Comparison  of  the  Coefficient  of  Variation  of  the  Ratio  of 
Ptheory/PACI  for  Concrete  Cylinder  Strength  Coefficients  of 
Variation  of  10%,  15%  and  20% 
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Table  5.11 


Comparison  of  the  Coefficient  of  Skewness  of  the  Ratio  of 
Ptheory/PACI  for  Concrete  Cylinder  Strength  Coefficients  of 
Variation  of  10%,  15%  and  20% 
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variation . 

The  mean  value  of  the  ratio  Ptheory/PACI  increased  in 
the  compression  failure  region  of  the  interaction  diagram 
with  increasing  cylinder  strength  coefficient  of  variation. 
This  increase  may  be  explained  by  the  increased  mean 
concrete  strength  required  by  ACI  318-71  Section  4.2  to 
account  for  the  increased  coefficient  of  variation.  There 
was  no  significant  increase  in  the  theoretical  strength  in 
the  tension  region  due  to  the  increased  mean  concrete 
strength.  Again  this  may  be  explained  by  the  compression 
failures  depending  on  the  concrete  strength  and  the  tension 
failures  depending  on  the  steel  strength. 

As  the  coefficient  of  variation  of  the  concrete 
cylinder  strength  was  increased  the  overall  coefficient  of 
variation  of  the  ratio  Ptheory/PACI  increased.  Again  the 
increase  in  overall  coefficient  of  variation  was  greater  in 
the  compression  region  of  the  interaction  diagram  where  the 
concrete  strength  has  more  effect  on  the  cross  section 
strength  than  in  the  tension  failure  region. 


5^3_j_4  The  Effect  of  The  Variables  Studied 

The  effect  of  the  variation  in  concrete  strength,  steel 
strength,  cross  section  dimensions  and  location  of  the 
reinforcing  steel  was  determined  for  the  12  in.  by  12  in. 
cross  section  with  a  nominal  steel  percentage  of  1%.  A 
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coefficient  of  variation  of  concrete  cylinder  strength  of 
15%  was  used  for  this  study.  Each  variables  effect  was 
studied  with  all  other  variables  at  their  mean  value. 

Figure  5. 1C  gives  a  graphical  representation  of  the 
overall  variation  in  cross  section  strength  for  various  e/h 
values  for  each  of  the  variables  and  for  all  the  variables 
combined.  The  plot  of  standard  deviation  squared  vs.  e/h 
indicates  the  major  component  causing  variation  in  cross 
section  strength  in  the  compression  region  of  the 
interaction  diagram  is  the  variability  in  the  concrete 
strength.  The  effect  of  the  variability  in  the  concrete 
strength  becomes  minimal  in  the  tension  failure  region.  This 
may  be  explained  by  the  fact  that  the  full  strength  of  the 
concrete  is  not  utilized  in  the  tension  failure  region  such 
that  the  high  concrete  strengths  have  no  effect  on  the 
variability  of  the  cross  section  strength. 

The  effect  of  the  variability  in  the  steel  strength  on 
the  overall  cross  section  strength  variability  is  greater  in 
the  tension  region  where  the  steel  strength  controls  the 
cross  section  capacity.  The  effect  of  the  steel  strength 
variability  in  the  compression  failure  region  is  again 
minimal  due  to  the  concrete  strength  being  the  controling 
factor . 


The  effect  of  the  variability  of  the  concrete  strength 
and  the  steel  strength  are  about  the  same  at  the  balance 
point.  This  is  to  be  expected  since  the  failures  in  a 
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randomly  selected  sample  would  depend  on  the  concrete  and 
steel  strength  equally  at  the  balance  point. 

The  effect  of  the  variability  in  the  cross  section 
dimensions  and  the  location  of  the  steel  was  very  small  for 
both  compression  and  tension  failures.  The  most  significant 
effect  occured  for  the  cases  of  pure  axial  load  and  pure 
moment . 

The  total  variability  in  the  cross  section  strength  may 
be  closely  approximated  by  the  expression: 

V  =  V  Z  +  V  +  V  /  (5.1) 

t  tc  ts  td 

where  may  be  the  total  standard  deviation  or  coefficient 

of  variation  and  V  ,  Vto  and  Vtd  are  the  standard  deviation 
or  coefficient  of  variation  of  the  cross  sectional  strength 
if  only  the  concrete  strength,  steel  strength  or  the 
dimensions  are  varied  separately. 

5 . 4  Cross  Section  Strength 

Figures  5.11  and  5.12  are  plots  of  the  interaction 
curves  for  the  12  in.  and  24  in.  columns  based  on  a  modified 
log-normal  distribution  of  steel  strength.  The  mean  strength 
indicated  is  the  mean  strength  calculated  from  a  sample  size 
of  500  using  the  Monte  Carlo  Technique  and  the  theoretical 
calculation  of  cross  section  strength.  The  maximum  and 
minimum  strength  curves  are  also  calculated  from  the  Monte 
Carlo  calculations.  The  ACI  ultimate  strength  is  the  cross 
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Figure  5.11  Dispersion  of  Strengths  of  an  Eccentrically 

Loaded  12  in.  Square  Column 
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Figure  5.12 


Dispersion  of 
Loaded  24  in. 


Strengths  of  an  Eccentrically 
Sguare  Column 
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section  capacity  calculated  using  the  ACI  318-71  Building 
Code.  The  ACI  factored  strength  is  the  ACI  ultimate  strength 
divided  by  1.4/0. 7  corresponding  to  the  lowest  possible  load 
factor.  The  ACI  factored  strength  corresponds  to  the  normal 
service  load  conditions.  The  discrepancy  in  the  mean 
strength  of  the  24  in.  column  immediately  above  the  balance 
point  is  due  to  the  reinforcing  steel  placement  in  the 
column  cross  section.  As  a  result  of  the  reinforcing  steel 
at  the  centre  of  the  cross  section  shifting  from  compression 
to  tension  steel  the  capacity  of  the  section  appears  to 
increase  to  a  second  balance  point  but  the  first  downturn  of 
the  curve  is  not  a  true  representation  of  the  capacity  of 
the  cross  section. 

The  dispersion  of  cross  section  strength  is  plotted  at 
selected  values  of  e/h.  In  each  case  the  dispersion  of  the 
cross  section  strength  is  a  normal  distribution  for  the 
compression  failure  region  and  a  log-normal  distribution  for 
the  tension  failure  region.  At  the  balance  point  the 
dispersion  of  cross  section  strength  may  be  represented 
equally  well  with  either  a  normal  or  log=normal 
distribution . 

Table  5.12  is  a  comparison  of  the  mean  value  and 
coefficient  of  variation  of  Ptheory/PACI  for  a  sample  size 
of  5C0  for  the  cross  sections  of  12  in.  by  12  in.  with  1% 
steel  and  24  in.  by  24  in.  with  3%  steel. 


In  the  compression  failure  region  the  variation  in  the 
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strength  of  the  12  in.  and  the  24  in.  columns  are  similar 
but  the  mean  value  of  Ptheory/PACI  is  larger  for  the  12  in. 
column.  This  is  due  to  the  higher  dependence  on  the  concrete 
strength.  Since  the  ratio  of  the  mean  concrete  strength  to 
the  nominal  concrete  strength  is  higher  than  the  ratio  of 
the  mean  steel  strength  to  the  nominal  steel  strength  the 
capacity  in  the  compression  region  increases  for  decreasing 
steel  percentages  relative  to  the  ACI  capacity. 

The  variability  of  the  theoretical  strength  of  the  24 
in.  by  24  in.  column  is  greater  than  that  of  the  12  in.  by 
12  in.  column  in  the  tension  region.  This  may  be  due  to  the 
increase  in  steel  percentage  from  1%  to  3%.  Also  the  mean 
value  of  the  ratio  Ptheory/PACI  is  larger  in  the  24  in. 
column  in  the  tension  region  due  to  the  increased  steel 
percentage. 

Figures  5.13  through  5.16  are  cumulative  frequency 
plots  of  the  ratio  Ptheory/PACI  for  the  12  in.  and  24  in. 
columns  at  selected  e/h  values.  A  comparison  of  these  and 
similar  cumulative  frequency  plots  for  a  normal  and  a  log¬ 
normal  dispersion  of  the  ratio  Ptheory/PACI  and  the  data  in 
Tables  5.1  to  5.4  and  5.9  to  5.11  suggest  that  for 
compression  failures  the  dispersion  may  be  represented  by 
the  normal  distribution  and  for  tension  failures  the 
dispersion  may  be  represented  with  a  log-normal 


distribution . 
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Figure  5.13  Normal  Cumulative  Frequency  Plot  of  the  Fatio 

Ptheory/PACI  for  the  12  in.  Column,  e/h  =  C.10 
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Figure  5.14  Log-normal  Cumulative  Frequency  Plot  of  the 

Fatio  Ptheory/PACI  for  the  12  in.  Column, 

e/h  =  0.10 
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Figure  5.15  Normal  Cumulative  Frequency  Plot  of  the  Ratio 

Ptheory/PACI  for  the  24  in.  Column,  Pure  Moment 
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Table  5.12 


Comparison  of  the  Mean  Value  of  the  Ratio  Ptheory/PACI  for 
the  12  in.  and  24  in.  Columns 
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5 . 5  Calculation  of  Factors 

5*5.1  Based  on  1_  in  J.00  Under  strength 

Tables  5.13  and  5.14  are  tables  of  the  calculated  <}> 
factor  based  on  a  probability  of  understrength  of  1  in  100 
and  a  normal  dispersion  of  cross  section  strength  in  the 
compression  failure  region  and  a  log-normal  dispersion  of 
cross  section  strength  in  the  tension  failure  region.  Figure 
5.17  is  a  plot  of  the  $  factor  for  the  12  in.  and  the  24  in. 
columns  vs.  e/h  based  on  a  probability  of  understrength  of  1 
in  100.  The  $  factors  in  the  ACI  Code  are  related  to  a 
probability  of  understrength  of  1  in  10040. 
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Table  5.13 


The  Understrength  Factor  for  the  12  in.  by  12  in. 
Based  on  a  Probability  of  Understrength  of  1  in  100 


r 


T 


l 


e/h 


Mean  Value 


Std.  Dev. 


$  Factor 


0.0 

\ 

1 

« 

1.22497 

0.05 

1 

1 

1 

1. 17698 

0.10 

1 

1 

I 

1. 15575 

0.  15 

I 

1 

« 

1. 13985 

0.20 

1 

1 

I 

1.12937 

0.30 

I 

1 

a 

1 . 12239 

0.40 

1 

1 

1 .12039 

0.50 

! 

1 

1 

1.06030 

0.60 

I 

1 

i 

1.04139 

0.70 

1 

1 

1 

1 . 03100 

0.80 

i 

1 

a 

1.02598 

0.90 

s 

I 

1.02439 

1.00 

! 

8 

i 

! 

1 

1.02498 

1.50 

1.C3976 

0^3 

i 

f 

■ 

1 .01598 

i _ — «- 


0.18915 

\  0.78 

0.17842 

|  0.76 

I 

0.17108 

1  0.76 

0.16583 

\  0.75 

0.16208 

8  0.75 

0.15369 

|  0.76 

0.13342 

|  0.81 

0.10804 

1  0.83 

0.10107 

f  0.83 

0.10035 

J  0.8  2 

0.10187 

|  0.81 

0.10374 

I  0.81 

0. 10522 

o 

00 

• 

o 

0.10649 

|  0.82 

l 

1  0.79 

0.10759 

j_ _ — _ 

Column 


112 


Table  5.14 


The  Understrength  Factor  for  the  24  in.  by  24  in.  Column 
Based  on  a  Probability  of  Under strength  of  1  in  100 
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e/h 


Fiaure  5  17  The  understrength  Factor  4>  vs.  e/h  Based  on  a 

Probability  of  Understrength  of  1  in  100  for 
the  12  in.  and  24  in.  Columns 
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5.5.2  Based  on  Cornell-Lind  Procedure 


It  has  been  proposed  that  future  code  revisions  have 


factors  based  on  the  equation: 


-8aV  RaV 

R  yDe  R  >  Uy  eP  u 
R  —  u 


(5.2) 


or: 


<j>  R  >  XU 


where : 


a  “3aV 

$  =  YRe  R 


=  the  understrength  factor 


>  _  e«v 

A  =  Yue  u 


-  the  overload  factor 


yR  =  Ptheory/PACI 


8  =  safety  index 


=  3.5  for  probability  of  failure  of  1.1  x  10  ~4 


=  4.0  for  probability  of  failure  of  3.2  x  10 


-s 


V  =  the  variability  of  the  strength  or  resistance 

B 


Vu =  the  variability  of  the  loads 


a  =  0.75 
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The  a  value  is  used  to  allow  the  separation  of  the 
effects  of  the  variability  of  the  member  strength  and  the 
variability  of  the  member  loading. 

Tables  5.15  and  5.16  are  tables  of  the  <p  factor  for  the 
12  in.  and  24  in.  columns  based  on  the  above  equation.  In 
these  tables  the  $  factors  are  based  on  values  of  6  of  4.0 
for  compression  failures  and  3.5  for  tension  failures.  The 
lower  probability  was  used  for  the  compression  failures  due 
to  the  sudden  brittle  mode  of  failure. 

The  12  in.  square  column  cross  section  was  chosen  to 
display  a  large  variability.  For  this  column,  <t>  was 
G.78±0.03  throughout  the  entire  range  of  e/h  studied.  The  24 
in.  column  was  chosen  to  have  a  low  variability.  For  this 
column,  was  0.79±0.09.  The  large  variability  in  the  24  in. 
column  was  due  to  the  discrepancy  in  the  theoretical 
strength  discussed  in  Section  5.4. 
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Table  5.15 


The  Understrength  Factor  for  the  12  in.  by  12  in.  Column 

„  ,  ,  -3aV„ 
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Table  5. 16 


The  Understrength  Factor  for  the  24  in.  by  24  in.  Column 
Based  on  <j>  =  y  e  -|3otV’ 
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CHAPTER  VI 


SUMMARY  AND  CONCLUSIONS 

In  this  study  probability  models  were  developed  to 
describe  the  variability  of  the  major  variables  affecting 
the  strength  of  a  reinforced  concrete  section.  Based  on  data 
from  a  literature  search  the  concrete  strength,  cross 
sectional  dimensions  and  location  of  reinforcing  steel  were 
described  with  a  normal  distribution  as  described  in  Chapter 
IV.  The  steel  was  described  with  a  normal  and  modified  log~ 
normal  distribution  of  yield  strength  as  discussed  in 
Appendix  A. 

A  Monte  Carlo  study  was  performed  using  the  probability 
models  developed  to  determine  the  variability  in  the  cross 
sectional  strength  of  a  12  in.  square  and  a  24  in.  square 
tied  reinforced  concrete  column.  The  results  of  this  study 
show  that  the  variability  of  the  concrete  strength  is  the 
major  contributing  factor  to  cross  sectional  strength 
variability  in  the  compression  failure  region  and  the 
variability  in  the  steel  strength  is  the  major  contributing 
factor  to  cross  sectional  strength  variability  in  the 
tension  failure  region.  The  effect  on  the  overall  strength 
variability  of  the  dimensional  variability  and  the 
variability  in  the  location  of  the  steel  strength  were  found 
to  be  minor  compared  to  the  effects  of  the  concrete  and 
steel  strength  variability.  The  type  of  distribution  assumed 
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for  the  steel  strength  variability  was  found  to 
significantly  affect  the  overall  strength  variability  in  the 
tension  failure  region  only. 

The  $  or  understrength  factors  were  calculated  based  on 
a  probability  of  understrength  of  1  in  100  and  based  on  the 
first  order  second  moment  procedure  developed  by  Cornell  and 
Lind.  The  calculated  values  of  <j>  were  in  close  agreement 
with  those  used  in  the  ACI  318~71  Code  for  column  cross 
sections  but  significantly  different  for  the  case  of  pure 
bending.  This  suggests  that  the  <p  factors  used  in  the  ACI 
Code  are  adequate  and  may  be  conservative  for  rectangular 
tied  column  cross  sections  but  seem  to  be  unconservative  for 
bending  tension  failures. 
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APPENDIX  A 


V API ABILITY  IN  REINFORCING  STEEL 


Introduction 

The  three  main  sources  of  variation  in  steel  strength 

are: 

(1)  variation  in  the  strength  of  material, 

(2)  variation  in  area  of  the  cross-section 
of  the  bar,  and 

(3)  variation  in  the  rate  of  loading. 

The  variability  of  yield  strength  depends  on  the  source 
and  the  nature  of  the  population.  The  variation  in  strength 
within  a  single  bar  is  relatively  small,  while  the  in-batch 
variations  are  slightly  larger.  However,  variability  Of 
samples  derived  from  different  batches  and  sources  may  be 
high.  This  is  expected  since  rolling  practices  and  quality 
measures  vary  for  different  countries,  different 
manufacturers  and  different  bar  sizes.  Furthermore,  the 
cross-sectional  areas  vary  due  to  differences  in  the  setting 
of  the  rolls,  and  this  adds  to  the  variation.  Mill  tests  are 
generally  carried  out  at  a  rapid  rate  of  loading  (ASTM 
corresponds  to  1040  micro-in/in/sec)  and  have  the  tendency 
of  reporting  the  unstable  high  yield  point  rather  than  the 
stable  low  yield  point.  Since  the  strains  in  the  structure 
are  induced  at  a  much  lower  rate  than  the  mill  tests,  mill 
tests  tend  to  overestimate  the  strength  of  reinforcement. 
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hence  another  source  of  variation. 

An  examination  of  the  test  data  revealed  that  the  bars 
of  large  diameter  tended  to  develop  less  strength  (4,  24, 
55)  than  #3  to  #11  bars.  Thus,  for  the  purpose  of 
statistical  evaluation,  the  #14  and  #18  bars  were  studied 
separately  from  the  other  sizes.  Also  the  #2  bars  were  not 
included  in  this  study  because  of  their  rare  use  for 
structural  concrete. 

In  this  study  the  terms  Grade  40,  Grade  50  and  Grade  60 
refer  to  reinforcing  bars  with  minimum  specified  yield 
strength  of  40,  50  and  60  ksi,  respectively,  even  though  the 
bars  in  question  may  not  have  been  produced  according  to 
ASTM  or  CS A  specifications.  Only  data  for  deformed  bars  has 
been  included.  In  some  cases  data  for  cold-worked  bars  has 
been  considered  but  most  of  the  data  is  for  hot-rolled  bars. 

Variation  in  Steel  Strength 

Different  values  for  the  yield  strength  of  steel  may  be 
obtained  depending  on  how  it  is  defined.  The  static  yield 
strength  based  on  nominal  area  seems  to  be  desireable 
because  the  strain  rate  is  similar  to  what  is  expected  in  a 
structure  and  designers  use  the  nominal  areas  in  their 
calculations.  Most  mill  tests,  however,  are  conducted  with  a 
rapid  rate  of  loading,  and  the  strength  is  generally 
referred  to  actual  areas*  For  these  reasons  the  yield 
strength  corresponding  to  rapid  strain  rate  and  measured 
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area  is  discussed  in  this  section,  the  effects  on  this 
strength  of  variations  in  cross-sectional  area  and  rate  of 
loading  are  dealt  with  in  the  succeeding  sections. 

A  review  of  literature  on  steel  strength  showed  that 
the  coefficient  of  variation  was  in  general  in  the  order  of 
1%  to  4%  for  individual  bar  sizes  and  4%  to  7%  overall  for 
data  derived  from  any  one  source.  When  data  was  taken  from 
many  sources  the  coefficient  of  variation  increased  to  5%  to 
8%  for  individual  sizes  and  10%  to  12%  overall.  A  summary  of 
selected  studies  from  literature  (4,8,9,33,43)  is  shown  in 
Table  A-1. 

The  data  reported  by  Allen4  and  Julian33  on  Grade  40 
and  60  steel  bars  showed  close  agreement  with  a  normal 
distribution  (with  respective  mean  and  standard  deviation) 
in  the  range  approximately  5  to  95  percentile  but  differ 
from  the  normal  distribution  outside  this  range.  Some 
authors  have  suggested  other  types  of  distributions  such  as 
skewed  distribution  (9,54),  truncated  normal  (31)  and  Beta 
distribution  (20) .  These  suggestions  were,  however,  based  on 
a  particular  set  of  data  and  only  approximated  the 
distribution  of  the  population  from  which  the  data  was 
drawn.  Nonetheless,  they  suggest  that  the  yield  strength  is 
a  phenomenon  that  can  be  described  by  a  particular 
theortical  distribution  with  certain  limitations.  The  normal 
distribution  seems  to  correlate  very  well  in  the  vicinity  of 
the  mean  values  for  different  populations  of  yield  strength/ 
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Summary  of  Selected  Studies  on  Steel  Strength 
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but  it  is  a  crude  approximation  at  low  and  high  levels  of 
probability  where  the  steel  strength  distribution  curves 
tend  to  have  certain  minimum  and  maximum  values  rather  than 
following  the  theoretical  tails.  This  is  expected  since 
there  are  always  some  quality  controls  that  are  used  to 
attain  a  certain  minimum  yield  strength  with  the  result  that 
the  manufacturing  of  steel  is  not  truly  a  random  process. 
Furthermore,  certain  data  indicated  a  positive  skewness, 
particularly  when  derived  from  different  sources  and  mixed 
together.  Theoretically  a  log-normal  distribution  should  be 
a  better  fit  for  this  case  than  a  normal  distribution  since 
it  takes  into  account  the  skew  nature  of  the  data.  However 
the  logarithmically  distributed  values  of  yield  strength  at 
low  and  high  levels  of  probability  did  not  show  a 
significant  improvement  over  normally  distributed  values  of 
available  data.  Therefore,  it  was  decided  to  empirically 
establish  a  "modified"  log-normal  distribution  that  would 
yield  correlate  with  the  North  American  data  on  yield 
strength. 


Values  of  (fy»34ksi)  are  plotted  on  log-normal 
probability  paper  in  Figure  A-1  for  the  data  from  Julian33 
and  Allen4  for  Grade  40  reinforcing  bars  grouped  together. 
The  values  were  found  to  be  in  good  agreement  with  a  log“ 
normal  distribution  in  the  range  from  the  0.01  percentile  to 
the  99th  percentile.  The  modification  constant  of  34  ksi  was 
established  by  trial  and  error.  The  corresponding  frequency 


curve. 


of  the  grouped  data  and  the 


the  histogram 


CUMULATIVE  FREQUENCY 
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(fy-34)  ( ksi) 


Figure  A~ 1 


Steel  strength  Distribution  for  Grade  40 
Reinforcing  Bars 
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corresponding  normal  frequency  distribution  curve  are  shown 
in  Figure  A-2  for  the  purpose  of  comparison.  The  mean  value 
of  the  data  was  found  to  be  48.8  ksi  with  a  maximum  value  of 
66  ksi  and  a  coefficient  of  variation  of  10.7%. 

Similarly,  values  of  (£y»55ksi)  for  Grade  60 
reinforcing  bars  from  mill  tests  reported  by  Allen4  were 
found  to  be  log-normally  distributed  in  the  range  from  the 
0.01  percentile  to  the  99th  percentile  as  shown  in  Figure  A- 
3.  The  frequency  curves  and  histogram  for  Grade  60  steel  are 
shown  in  Figure  A-4.  The  mean  value  for  the  data  was  71.5 
ksi  with  a  maximum  value  of  90  ksi  and  a  coefficient  of 
variation  of  7.7%. 

In  both  cases  the  modified  log-normal  curve  is  a  better 
approximation  at  the  lower  end  of  the  curve  while  the  normal 
curve  is  better  at  the  high  end  of  the  curve. 

Variation  in  Steel  Cross-Sectional  Area 

The  actual  areas  of  reinforcing  bars  tends  to  deviate 
from  the  nominal  areas  due  to  the  rolling  process.  The 
designers  do  not  have  this  information  readily  available  to 
them,  and  hence  use  the  nominal  areas  in  their  calculations. 
For  this  reason  this  variation  should  be  incorporated  in  the 
strength  of  steel. 

variation  in  the  ratio  of  measured  to  nominal 
areas  (A  /AR)  was  studied  as  a  measure  of  the  variation  in 


_ 
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C  *  r  1  fX-X 

PDF  =  - t=  exp  [-  2 


yax/2TT 


)  ] 


X 


X  =  1.19456,  a  =  0.14112 

X 

c  =  0.43429 

y  =  (fy  -  34  ksi)  ,  x  =  LOG^y 


NORMAL  DISTRIBUTION 

X  =  48.8  ksi,  a  =  5.506  ksi 

X  =  fy  ksi 


Figure  A-2 


Probability  Density  Function  for  Grade  40  Bars 
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(fy- 55)  (ksi) 


Figure  A~3 


Steel  Strength  Distribution  for  Grade  60 
Feinforcing  Bars 
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Figure  A-4  Probability  Density  Function  for  Grade  60  Ears 


the  cross-sectional  area  of  reinforcing  bars.  The  values  of 
^e/An  are  reproduced  from  available  literature  (4,  8,  43)  in 
Table  A°1.  Table  A- 1  indicates  that  the  data  reported  by 
Baker8  for  Grade  60  steel  demonstrates  high  mean  value  and 
coefficient  of  variation.  Such  values  cannot  be  explained  in 
definite  terms.  It  is  possible  that  the  collected  data 
contained  a  good  percentage  of  values  from  mills  with  old 
rolls  that  increased  the  mean  and  coefficient  of  variation. 
Furthermore,  British  rolling  practice  may  differ  from 
Canadian  practice.  For  these  reasons,  these  values  were  not 
included  in  the  analysis. 

The  ratios  of  Ae/An  from  tests  on  Grade  40  and  60 
reinforcing  bars,  manufactured  in  Canada  (Study  No.  1  and  3 
in  Table  A-1),  were  plotted  on  normal  probability  paper. 
These  values  exhibited  close  agreement  in  the  range  from  the 
5th  to  the  95th  percentile  for  Grade  40  steel  and  from  the 
2nd  to  the  98th  percentile  for  Grade  60  bars  with  a  normal 
distribution.  When  the  values  for  both  studies  were  combined 
they  resulted  in  a  normal  distribution  in  the  range  between 
the  4th  and  99th  percentile  with  a  mean  value  of  0.988  and  a 
coefficient  of  variation  of  2.4%.  The  effect  of  such  a  small 
coefficient  of  variation  is  not  large  enough  to  have  any 
significant  effect  on  the  coefficient  of  variation  of  Asfy. 
For  this  reason  a  single  value  for  Ae/An  seems  to  be  more 
appropriate.  Allen4  has  suggested  a  value  of  0.97  for  A0/An . 
This  seems  to  be  a  conservative  estimate  of  the  average 
values  of  Ae/An  shown  in  Table  A*-1r 


and  close  to  ASTM 


.. 
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rolling  tolerances  that  allow  an  average  ratio  as  low  as 
0.965  and  a  minimum  single  value  up  to  0.940. 

Effect  of  Rate  of  Loading 

The  apparent  yield  strength  of  a  test  specimen 
increases  as  the  strain  rate  or  the  rate  of  loading 
increases.  Since  mill  tests  on  steel  specimens  are  generally 
carried  out  at  much  greater  strain  rates  (approximately  1040 
micro- in/in/sec)  than  encountered  in  a  structure,  they  tend 
to  overestimate  the  yield  strength.  A  strain  rate  of  1 
in/in/sec  may  increase  the  yield  strength  of  Grade  40  steel 
as  much  as  50%  over  the  static  yield  strength  (34) . 

Tests  conducted  on  steel  coupons  of  A36,  A441  and  A514 
steel  (51)  demonstrated  a  yield  strength  reduction  more  or 
less  of  the  same  value  for  all  types  of  steel  with  decrease 
in  the  rate  of  strain.  The  equation  developed  by  Fao5*  on 
the  basis  of  these  tests  gives  values  of  static  yield 

strength  that  are  4.8  ksi  and  3.4  ksi  less  than  the  yield 

strengths  obtained  at  cross-head  speed  of  1000  and  200 
micro~in/in/sec  respectively.  NRC  tests  on  Grade  40  bars  (4) 
showed  a  reduction  of  approximately  3  ksi  in  the  mean  yield 

strength  when  speed  of  the  testing  machine  was  dropped  from 

208  m i cr o— i n/i n/sec  to  static.  This  value  correlates  well 
with  the  one  obtained  from  Bao's  equation.  Similarity,  for 
Grade  40  bars,  it  has  been  shown  at  the  University  of 
Illinois  (34)  that  the  difference  between  the  yield  strength 
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at  a  strain  rate  of  1040  micro- in/in/sec  and  the  strength  at 
a  strain  rate  of  20  micro-in/in/sec  is  about  9%  or  4  ksi. 
ETH  tests  (36)  for  high  strength  reinforcement  demonstrated 
a  reduction  of  3  ksi  for  static  conditions. 

For  evaluation  of  the  static  yield  strength  from  mill 
tests,  Allen4  has  suggested  a  decrease  of  4  ksi.  This  value 
seems  to  be  a  reasonable  estimate  for  the  available  test 
data . 

Effect  of  Bar  Diameter 

The  strength  of  steel  tends  to  vary  across  the  cross- 
section  of  a  reinforcing  bar  with  the  highest  strength  near 
the  outside  of  the  bar.  This  is  probably  due  to  the  cold¬ 
working  of  circumferential  sections  of  bars  during  the 
rolling  process.  Thus  the  mean  yield  strength  is  expected  to 
decrease  with  increasing  diameter.  The  variation  of  the  mean 
yield  strength  with  size  is  plotted  in  Figures  A-5  and  A-6. 
The  data  shown  in  the  figures  were  taken  from  several  test 
series  for  Grade  40  and  Grade  60  reinforcement 
(4,44,8,9,24).  For  bars  with  relatively  small  diameter  the 
effect  of  this  variation  is  small  and  not  clearly 
established.  For  large  diameter  bars  such  as  #14  and  #18 
this  effect  becomes  more  significant.  In  addition,  the  A STM 
specifications  allow  the  use  of  small  specimens  machined 
from  samples  of  large  diameter  bars  for  testing  purposes.  A 
specimen  machined  to  a  smaller  diameter  from  a  guarter- piece 
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of  a  full  size  bar  tends  to  show  higher  yield  strength  than 
the  bar  itself  (24) .  Since  some  manufacturers  may  use  these 
tests  as  a  measure  of  quality  control,  the  #14  and  #18  bars 
tend  towards  a  higher  probability  of  passing  through  quality 
controls  without  developing  the  required  strength. 

An  extremely  limited  amount  of  data  is  available  for 
#14  and  #18  bars.  Tests  on  Grade  40,  #14  bars  carried  out  by 
Allen4  showed  that  the  mean  yield  strength  of  #14  bars  was 
44  ksi,  a  15%  decrease  from  the  strength  of  #3  to  #11  bars 
produced  by  the  same  manufacturer.  Some  data  has  been 
reported  by  Gamble24  for  #14  and  #18  bars  of  Grade  60  steel. 
The  mean  yield  strengths  were  60  ksi  for  #14  and  55  ksi  for 
#18  bars.  These  strengths  were  referred  to  the  nominal 
areas.  Using  the  mean  yield  strength  of  Grade  60,  #3  to  #11 
bars  as  71.5  ksi  (as  per  Study  No.  3  in  Table  A~1)  and  a  3% 
adjustment  for  the  deviation  from  the  nominal  area,  the 
reduction  in  strength  is  approximately  13%  for  #14  bars  and 
21%  for  #18  bars.  This  comparison  is,  however,  not  truly 
justified  since  the  data  for  both  studies  was  not  drawn  from 
the  same  source.  Nonetheless,  it  strongly  indicates  the 
understrength  of  #14  and  #18  bars.  Until  more  data  is 
available,  it  seems  reasonable  that  the  yield  strength  of 
#14  and  #18  bars  should  be  reduced  at  least  15%  below  the 
yigj_d  strength  of  reinforcing  bars  with  smaller  diameter. 
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Summary 

The  modified  log-normal  distribution  curves  shown  in 
Figures  A- 1  through  A-4  seem  to  correlate  well,  particularly 
near  the  lower  tails  of  the  curves,  with  the  available  North 
American  data  for  Grade  40  and  Grade  60  reinforcing  bars. 
The  Probability  Density  Function  for  these  curves  can  be 
calculated  using  the  following  equation: 

PDF  =  c 

y a  »/  2n 

x 

where : 


c  =  0.43429 


y 

y 


34  ksi  for  Grade  40  bars 

55  ksi  for  Grade  60  bars 


x  = 


Log 

10 


y 


x  =  1.14482  for  Grade  40  bars 


x  =  1.19456  for  Grade  60  bars 

a  =  0.14866  for  Grade  40  bars 
x 

a  =  0.14112  for  Grade  60  bars 

x 

The  mean  yield  strength  of  the  selected  data  was  found 
to  be  48.8  ksi  (c.o.v.  =  10.7%)  for  Grade  40  bars  and  71.5 
ksi  (c.o.v.  =  7.7%)  for  Grade  60  bars.  The  modification 

constants  were  empirically  established  and  found  to  be  34 
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ksi  and  55  ksi  for  Grade  4C  and  Grade  60  steel  respectively. 

A  value  of  0.97  for  the  ratio  Ae/AR  seems  to  be 
reasonable  to  account  for  deviations  from  the  nominal  areas. 
Similarly,  for  the  evaluation  of  the  static  yield  strength, 
at  least  4  ksi  should  be  deducted  from  the  yield  strength- 
obtained  in  mill  tests  or  at  high  strain  rates  allowed  by 
ASTM  specifications. 

When  calculating  the  yield  strength  of  #14  and  #18 
reinforcing  bars  from  the  strength  of  bars  of  smaller  sizes 
at  least  a  15%  reduction  should  be  used  to  account  for  the 
effect  of  the  large  diameter. 


APPENDIX  B 


COLUMNS  STUDIED 

This  appendix  contains  the  details  of  the  two  major 
columns  studied.  Tables  B- 1  and  B-2  are  tables  of  the 
properties  of  the  12  in.  and  24  in.  columns  respectively. 
Figures  B-1  and  B-2  are  diagrams  of  each  column  showing  the 
designer's  properties  and  the  mean  values  of  the  column 
properties  used  in  the  Monte  Carlo  calculations. 
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Table  B- 1 


Properties  of  the  12  in.  Column  Assumed  in  the  Calculations 

Specified  Mean  In-  o  C.C.V. 

situ 


Material  Strengths 


Concrete  Strength 

3000 

psi. 

3712 

psi. 

0.17 

Steel  Yield  Strength 

40 

ksi . 

48.8 

ksi . 

1.41 

ksi . 

Dimensions 

b ,  h 

12.00 

in . 

12.06 

in . 

0.280 

in. 

d 

9.75 

in . 

9.51 

in . 

0.166 

in.  *=-- 

d' 

2.25 

in . 

2.55 

in. 

C.166 

in.  --- 

*s 

1.76 

sg .in . 

on  cm 

As 

0.88 

sg .in . 

0C~j  ar1  c 

s 

12.00 

in . 

12.00 

in . 

—  - 

<c  cs  e 

b",  d" 

9.00 

in . 

8.47 

in . 

0.166 

in.  — - 

A" 

0.1 1 

sg  .in . 

0.11 

sg .  in 

• 

Individual  Longitudinal  Steel  Bars 


ASB  (1) 

to  ASB  (4) 

0.44 

sg  .in . 

0.44 

sg.in.  --- 

DS  ( 1)  , 

DS  (2) 

2.25 

in . 

2.55 

in . 

0.166  in . 

DS(3)  » 

DS  (4) 

9.75 

in . 

9.51 

in . 

0.166  in. 

ta  lt. 
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Table  B- 2 

Properties  of  the  24  in.  Column  Assumed  in  the  Calculations 

Specified  Mean  In-  a  C.O.V. 

situ 

Material  Strenqhts 


Concrete  Strength 

30 00  psi. 

3712 

psi . 

C.  176 

Steel  Yield  Strength 

4  0  k  si . 

4 

CO 

• 

00 

ksi 

1, 

.41 

ksi . 

«T3  c 

Dimensions 

b,  h 

24 

o 

o 

• 

in . 

24 

.06 

in  . 

0. 

,280 

in . 

d 

21 

.30 

in . 

21 

.Cl 

in . 

0. 

.166 

in . 

d* 

2 

.70 

in . 

3 

.05 

in . 

0. 

,166 

in . 

mo  mo 

As 

18 

.72 

sq.in. 

- 

— 

3  CEE-  ru 

a,a6m 

7 

o 

CO 

• 

sg  .in . 

- 

— 

—  - 

S 

12 

o 

o 

• 

in . 

12 

o 

o 

in . 

b",  d" 

20 

.50 

in . 

19 

.87 

in . 

0. 

166 

in . 

«o 

A”  0.20  sq.in.  0.20  sq.in.  --«*  — - 

Individual  Longitudinal  Steel  Bars 
ASB(1)  to  ASB(12)  1.56  sq.in.  1.56  sq.in.  — 

DS ( 1 )  to  DS (5)  2.70  in.  3.05  in.  0.166  in.  --- 

DS (6)  ,  DS (7)  12.00  in.  12.07  in.  0.993  in. 

21.30  in.  21.01  in.  0.166  in.  — - 


DS (8)  to  DS (12) 


148 


A 

c 

ID 

v 


c 

C\l 


f  =  40,000  psi 

f  =  3,000  psi 

A$  —  2#  6  bars 
A;  -  2#  6  bars 
A"  -  #3  @12  in 


12  in.  x  12  in.  COLUMN 


f  =  40,000 

fy  =  40,000  psi 

f'  =  3,000  psi 

A  —  7#  1 1  bars 

S 

A'  —  5#  1 1  bars 

S 

A''  —  #4  @12  in. 
in  sets  of  3 


24  in.  x  24  in.  COLUMN 


Figure  B-1 


Nominal  or  Designer's  Properties  of  the  12  in. 
and  24  in  Columns 


- 
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fy  =  48,800  psi 

f  =  3712  psi 
A  —  2#  6  bars 

S 

A'  —  2#  6  bars 

S 

A"-  #3  @12  in 

S 


12  in.  x  12  in.  COLUMN 


fy  =  48,800  psi 

fc  =  3712  psi 

A  —  7#  1 1  bars 

S 

A;  -  5#  1 1  bars 

A"-  #4  @12  in. 
in  sets  of  3 


24  in.  x  24  in.  COLUMN 


Figure  B—  2  M©an  Valu©s  of  fh©  Piropoirtiss  of  fh©  12  in. 

24  in  Columns 


and 


APPENDIX  C 


FLOW  DIAGRAMS  OF  THE  MONTE  CARLO  PROGRAM 

This  appendix  contains  detailed  flow  diagrams  of  the 
complete  Monte  Carlo  Program  including: 

The  Main  Program 
Subroutine  PROP 
Subroutine  ACI 
Subroutine  ASTEEL 
Subroutine  CURVE 
Subroutine  THMEAN 
Subroutine  THEORY 
Subroutine  AXIAL 
Subroutine  F STEEL 
Subroutine  RANDOM 
Subroutine  STAT 
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MONTE  CARLO  PROGRAM 


, 


152 


£ 


II  1  MQ 

V 

JJ  1 ,  IMo 

I  =  1,  NV 

I 


4- 


SD 

=  STDV  (1) 

CONST  =  FCONST  (1) 

RM 

=  RMEAN  (1) 

ITP  =  ITYPE  (1) 

CALL  RANDOM 


FY  =  0.97  (X  (2)  -4000) 

5 

BB  =  X  (3)  Dll  =  X  (6) 

H  =  =  X  (4)  DC  =  X  (7) 

B11  =  X{5)  DD  =  X  (8) 

3Z 

I  =  1,  NB 


Nl  =  I  +  8 


DS  (1)  = 

=  X  (Nl) 

(!) 


1  55 


156 


SUBROUTINE  PROP 


157 


SUBROUTINE  ACI 


A 


' 


158 


I 


1  59 


£ 


P  (J)  =  FCCONC  +  FST 


COMPM  =  FCCONC*(DD  -  B1  ■  C/2 


SM  =  0.0 


I  =  1,  NB 


SM  =  SM  +  SBM  (I) 


BM  (J)  =  COMPM  +  SM  -  P  (J)*(DD  -  H/2) 


EOH  (J)  =BM  (J)/(P(J)*H)  4- 


YES 


YES 


WRITE  PO,  PB,  BMO 


© 


SBM  (I)  =  (FSS  (I)  -  FCS  (l)*ASB  (l))*DD  -  DS  (I)) 


WRITE  ACI  INTERACTION 
DIAGRAM 


160 


SUBROUTINE  ASTEEL 


161 


SUBROUTINE  CURVE 


DEFINE  EOH1  (1) 

1  =  1,  NN 

< 

r 

i  = 

1,  N 

YES 


CONTINUE 

r 

•i 

NP 

=  1 

1 

EOHB  = 
Ml  =  M 

EOH  (1) 

YES  ^ 

M  = 

NP— 2 

164 


§ 


© 


I  =  1,  JJ 

m 


P  (II)  =  P(ll)  +  COE  (I)* [EOH1  (ll)**(l  -  1)] 


CONTINUE 


NNN  =  II 


I  =  NP.  N 


DIFE  =  3.0  —  EOH  (I) 


CONTINUE 

n 

► 

1 

NM  =  N 

-  NP  +  1 

M  = 

Ml 

N  =  1  -  1 

M  =  NM  —  2 

165 


©  ® 


YES 


SUMIP  =  ANS  (4)  -  SUM 


SUMIP  >  0.0 
YES 


SUM  =  ANS  (4) 
COE  (1)  =  ANS  (1) 


J  = 

U 

t _ 

COE  (J  +  1)  =  B  (J) 

_ 

[ 

JJ  =  I  +  1 


T 


CONTINUE 


I 


II  =  NNN,  NN 


BM  (II)  =  0.0 

~T~ 


I  =  1,  JJ 
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168 


SUBROUTINE  THMEAN 


169 


SUBROUTINE  THEORY 


1  70 


172 


SUBROUTINE  AXIAL 


A 


B 


173 


© 


P11  =  (2*(B1 1  +  D1 1)*AS1 1  +  ASC*S)/(B11*D11*S) 


E50H  =  0.75*P11  VB11/S 

E50U  =  (3  +  0.002*FC)/  (FC  -  1000) 


z  = 
zz 

EU 

0.5/(E50H  +  E50U  -  EO 
=  0.5/(E50U  -  EO) 

=  2*(E50H  +  E50U)  -  EO 

DX  =  C/10 

1  =  1,10 

%  . 

i  A,=l  

* 

X  (1)  =  C  -  AI*Dx  +  Dx/2 

E  (1)  =  PHI*X(I)  +  ECO 

B  (1)  =  BB 

-® 


B  (I)  =  B11 


174 


175 


■ 


176 


177 


SUBROUTINE  F  STEEL 


178 


SUBROUTINE  RANDOM 


1  79 


SUBROUTINE  ST  AT 


18C 


C 


181 


183 


184 


APPENDIX  D 


LISTING  OF  THE  MONTE  CAPLO  PPOGPAM 

This  appendix  contains  a  complete  listing  of  the  Monte 
Carlo  Program.  The  modified  IBM  Subroutine  MULT  to  MULTR  is 
also  listed.  The  listing  includes: 

The  Main  Program 
Subroutine  THMEAN 


Subroutine 

A  Cl 

Subroutine 

ASTEEL 

Subroutine 

PROP 

Subroutine 

CURVE 

Subroutine 

RANDOM 

Subroutine 

THEORY 

Subroutine  AXIAL 
Subroutine  FSTEEL 
Subroutine  STAT 
Subroutine  TMULTR 


185 


186 


$  LI ST  F2 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 


ON  *  PEI  NT* 

£******************************«********************** 

COMMON  N.EOH1  (1 3)  , PO , B MO , DCS , DTS 

COMMON  FC,FY,ES,BB,H,DC,DD,AS,  AS  1 1 , B 1  1  ,D1  1  ,  S  ,  C  ,  27  ,  NED 
COMMON  PH  I , EO , J , 2 ,  ECC ,  EY  ,  FCONCC  (20)  ,  ASC,  EOH  (40)  ,  FCONST  (25) 
COMMON  X  (16  000)  ,  EC (20)  ,  B  (2  0)  ,  P  (4  0)  ,BMM  (1000)  ,  BM (4 0)  ,  FCS  (20) 
COMMON  FST ,  E  (20 )  ,NB,DS  (20)  ,ASB(2  0)  ,FS5  (2  0)  ,SBM  (20) 

DIMENSION  PACI  (15), PTH  (15, 2000), RP  (15,  200  0) 

DIMENSION  FMFAN  (25)  ,STDV  (2  5)  ,ITYPE  (25) 

DIMENSION  CLASS  ( 3  1 )  , CFREQ ( 3 1 )  , NG BOSS  (3 1 ) 

C  BEAD  QUANTITIES  NEEDED  FOE  MONTE  CARLO  SIMULATIONS 
READ  (5 , 500)  NV,NS,FY1 
500  FORMAT  (215, F9. 2) 

BEAD  (5,511)  RMEAN1,RMEAN2,IY,NRU 

511  FORMAT  (2F 1 5 . 5 , 21 1 0) 

C  READ  NOMINAL  PROPERTIES 

CALL  PROP  (NS) 

C  BEAD  STATISTICAL  PROPERTIES  OF  VARIABLES 
DO  5  1=1, NV 

READ  (5,510)  RMEAN(I)  ,STDV(I)  ,  FCONST  (I)  ,ITYPE(I) 

510  FOFMAT(3F15.5,I5) 

5  CONTINUE 

C  WRITE  STATISTICAL  PROPERTIES  OF  VARIABLES 
WRITE  (6,512) 

512  FORMAT  (’ 1* ,//////27X, •DISTRIBUTION  PROPERTIES  OF  VARIABLES') 
WRITE  (6,519)  NS , NRU 

519  FORMAT  (3 1 X , •  ( •  ,  14  ,  • S IM • , 1 5, • )  •//) 

WRITE  (6,513) 


513 

1000 

514 

515 

516 
522 
520 


MEAN-VALUE  ST D-D E VI  AT ION 


FORMAT  (16X, ' 

DO  1000  1=1, NV 

WRITE  (6,514)  RMEAN  (I)  ,STDV  (I)  , FCONST (I)  ,ITYPE  (I) 

FORMAT  (16X,3F15. 5,15) 

WRITE  (6,515) 

FORMAT  (////2G X , '  FC  (MEAN- VALUE)  FY  (MEAN-VALUE)  '/) 

WRITE  (6,516)  RMEAN 1 , RMEAN2 
FORMAT  (20X, 2F15. 5) 

WRITE  (6,522)  FY 1 
FORMAT  (///21X,'FY  LIMIT=* ,F9.2) 

WRITE  (6,520)  IY 
FORMAT  (///2 IX, 1 ISEED=' ,110) 

CALCULATE  THE  ACI  INTERACTION  DIAGRAM 
CALL  ACI  (NS) 

FIT  A  POLYNOMINAL  TO  THE  ACI  INTERACTION  DIAGRAM 
CALL  CURVE 

WRITE  THE  ACI  INTERACTION  DIAGRAM  AFTER  THE  CURVE  FIT 

(6, 100) 

(' 1 ’ ,//////23X, ' ****ACI  INTERACTION  DIAGRAM****') 
(6,519)  NS, NRU 
(6,517) 

(30X,*  (AFTER  CURVE  FIT)'//) 

(6,103) 

(19X,'P(J)  LBS'  ,6X,'M(J)  LB-IN'  ,7X,'E0H(J)  •) 

DO  6  0=1,13 

WRITE  (6,104)  P  (J)  ,BM  (J)  ,EOHl  (J) 

FORMAT  (/1 6X ,  3E1  5. 7) 

PACI  (1)=PO 


FCONSTANT  TYPE'/) 


100 

517 

103 
6 

104 


WRITE 

FORMAT 

WRITE 

WRITE 

FORMAT 

WRITE 

FORMAT 


56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

11C 

111 

112 

113 

114 

115 


187 


PACI  (15) -8R0 
DO  1  1=1,13 
1  PACI  (1  +  1)  =P  (I) 

C  CALCULATE  MEAN  THEORY  INTERACTION  DIAGRAM 
CALL  THMEAN  ( F.  MEAN  ,  P  M  E  A  N  1  ,  RME  A  N  2  ,  NS) 

C  FIT  A  POI.YNCMJNAL  TO  THE  THEORY  INTERACTION  DIAGRAM 
CALL  CURVE 

C  WRITE  MEAN  THEORY  INTERACTION  DIAGRAM  AFTER  CURVE  FIT 
WRITE  (6,101) 

101  FORMAT  (' 1' ,//////1QX, • ****MEAN  THEORY  INTERACTION  DIAGRAM****') 
WFITE  (6,519)  NS , NPU 
WRITE  (6,517) 

WRITE  (6,103) 

DO  3  J=1 , 13 

3  WRITE  (6,104)  P  (J)  ,BK(J)  ,EOH1  (J) 

PTK  (1  ,  1)  =  PO 
PTH  (15, 1)  =EMO 
DO  15  1=1, 13 
15  PTH  (1+1,1)  =P  (I) 

WFITE  (6,105) 

105  FORMAT  (//1 8X , • MEANT H/ACI '  , 2X  ,  *  EOH’ ) 

C  CALCULATE  AND  WRITE  RATIO  MEAN  THEORY/ACI 
DO  17  1=1,15 
IF  (I.EQ.1)  EOH 2=0 . 0 
IF  (I.GT.1)  EOH2=EOHl  (1-1) 

IF  (I.EQ.15)  EOH2=99 . 9  9 
RP  (I  ,  1 )  =FTH  (I,  1 )  /PACI  (I) 

17  WRITE  (6,106)  RP(I,1),EOH2 

106  FORMAT  ( 1 6 X, 2F 1 0 . 5) 

C  MONTE  CARLO  CALCULATION  OF  THEORETICAL  STRENGTH 
DO  4  J J= 1 , NS 
DO  10  1=1, MV 
SD  =  STDV  (I) 

RM=RME AN (I) 

CONST-FCONST  (I) 

ITP=ITYFE  (I) 

CALL  RANDOM  (I Y , S D , RM , CONST, ITP, V) 

X(I)=Y 
10  CONTINUE 

IF  (X(1)  .IE.  (RMEAN(1)-3.3*STDV  (1)))  X  ( 1 )  =RMEAN  ( 1)  -  3 . 3*  STD  V  ( 1 ) 
FC=X (1) 

IF  (X  (2)  .GT.FY1)  X  (2)  =FY  1 
FY=  (X(2)-4000.0)*0.97 
BB=X  (3) 

H=X  (4) 

B1 1  =  X  (5) 

Dll  =X  (6) 

DC=X  (7) 

DD  =  X  (8) 

DO  2  1=1, NB 
Nl=l+8 

2  DS (I ) =X  (NI) 

C  CALCULATE  THEORETICAL  INTERACTION  DIAGRAM 
CALL  THEORY 

C  FIT  A  POLYNOMIAL  TO  THE  THEORY  INTERACTION  DIAGRAM 
CALL  CURVE 
PTH  (1,  JJ)  =PO 
PTI!  (15, JJ)  =BMO 
DO  9  1=1,13 
9  FTH  (1+1, OJ) -P(I) 


116 

1 17 

1 18 

119 

120 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 

161 

162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 


188 


C  CALCULATE  RATIO  PTHEOPY/PACI  AT  SPECIFIED  E/H  VALUES 
DO  44  1=1  , 15 

R  P  (1  ,  J  J)  =PTH  (I,JJ)/PACI  (I) 

44  CONTINUE 
4  CONTINUE 

C  STATISTICAL  ANALYSIS  OF  PTHEOSY/PACI  FOR  EACH  SPECIFIED  E/H  VALUE 
DO  45  1=1,15 
DO  8  J J= 1 , NS 
L=J  J 

8  X  (I,)  =RP  (I ,  JJ) 

N  =  NS 

CALL  STAT  (R K, SD , COV, COS , COK , RMI N , RH A X , U N 2 , UH 3 , UM4 , I BA X, I HI N , RMED , 
1 CLASS , CFREQ ,NGROSS) 

IF  (I.EQ.1)  ECH2=0 . 0 
IF  (I.GT.1)  EOH2=EOH1 (1-1) 

IF  (I.EQ.15)  EOH2  =  99 . 9  9 
WRITE  (6,518) 

518  FORMAT  ( *  1 • ,////// 3 1 X, ******* +THEORY/A Cl ***♦***') 

WRITE  (6,22)  EOH2 

22  FORMAT  (,/3  1  X  ,*********  EOH=  ',  F5 . 2  ,'  ******** ) 

WRITE  (6,519)  NS , NRU 
WRITE  (6,16) 

16  FORMAT  (3 1 X,  • <*>  STATISTICAL  EVALUATION  <*>') 

WRITE  (6,20) 

20  FORMAT  (/I 6X ,' MEAN- VALUE  SD-DEVIATION  CO-VARIATION  CO-SKEWNESS  K 
1 URTOSIS' ) 

WRITE  (6,25)  EM, SD, COV, COS, COK 
25  FORMAT  (1 3 X , 4F1 3 . 5 , FI  0 . 5) 

WRITE  (6,30) 

30  FORMAT  (/ 1 7 X , ' MI N- V ALU E  (SIMULN  NO.)  MAX-VALUE  (SIMULN  NO.) 

1  MEDIAN ' ) 

WRITE  (6,35)  RMIN,IMIN,RMAX,IMAX ,RMED 

35  FORMAT  (1  3X,F1 3. 5, 1 13, F13. 5,  1 1 3,  F 10 . 5) 

WRITE  (6,36) 

36  FORMAT  (/25X, ' MOMENTS  ABOUT  THE  MEAN') 

WRITE  (6,37) 

37  FORMAT  (18X, '2ND- MO ME NT  3RD- MOMENT  4TH-MOMENT ' ) 

WRITE  (6,38)  UM2,UM3,UM4 

38  FORMAT  (16X,3E15.7) 

WRITE  (6,39) 

39  FORMAT  (23X, • CUMULATIVE  FREQUENCY  TABLE') 

WRITE  (6,40) 

40  FORMAT  ( 1 9 X , *  CL ASS- NO .  UPPER-LIMIT  %CUM-FREQ.  G ROSS- NUM BER • ) 

DO  45  111=1,31 

45  WRITE  (6,50)  III , CLASS  (III) , CFREQ (III)  , NGROSS  (III) 

50  FORMAT  (15X,I13,2F13.5,1I13) 

WRITE  (6,1505) 

1505  FORMAT  (' 1 ' ,//////26X, ' TOTAL  POPULATION:  PHI  FACTORS') 

WRITE  (6,519)  NS, NRU 

C  STATISTICAL  ANALYSIS  CF  PTHEORY/PACI  FOR  ALL  E/H  VALUES  COMBINED 
L=0 

DO  7  1=1,15 

IF  (I.EQ.1)  EOH2=0 . 0 

IF  (I.GT.1)  EOH2=EOH1 (1-1) 

IF  (I.EQ.15)  EOH2=99 .99 
WRITE  (6,1515)  EOH2 

1515  FORMAT  (/26X,'***EOH=' ,F5. 2, '***'/) 

WRITE  (6,  1500)  (RP(I,JJ),  JJ  =  1,NS) 

1500  FORMAT  (21 X, 5F10. 5) 

DO  7  JJ- 1 , NS 
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177 

178 

179 

180 

181 

182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 

201 

202 

203 

204 

205 

206 

207 

208 

209 

210 

211 

212 

213 

214 

215 

216 

217 

218 

219 

220 

221 

222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 


1  89 


L=L  +  1 

7  X (L) =  RP  (I, JJ) 

N  =  L 

CALL  STAT  ( R M , S D , COV , COS ,COK , RMI N, RMAX , 0 M2, UM3 , UM4 ,IMAX,IMIN,RMED, 
1CI,ASS,CFREQ,NGROSS) 

DC  11  1=1,15 
IHAX=IMAX-NS 
IF  (IMAX.LE.O)  GO  TO  12 

11  CONTINUE 

12  IMAX=IMAX+NS 
DO  13  1=1,15 
IMIN=I MIN-NS 

IF  (IMIN.LE.O)  GO  TO  14 

13  CONTINUE 

14  IMIN=IMIN+NS 
WRITE  (6,21) 

21  FORMAT  (' 1' ,//////26X, •  <*> 

(6,519)  NS , NRU 

(6,20) 


FORMAT 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

DO  55 

WRITE 


TOTAL  STATISTICAL  EVALUATION  <*>') 


(6,25) 

(6,30) 

(6.35) 

(6.36) 

(6.37) 

(6.38) 

(6.39) 

(6.40) 

1=1,31 

55  WRITE  (6,50) 

1600  CONTINUE 

WRITE  (6,1900)  JJ 

1900  FORMAT  {' 1  •  ,/20 X ,'***', 15 ,'***' ) 

STOP 

END 

C***************************************************** 


RM,SD,COV,COS,COK 

RMIN,IMIN,RMAX,IMAX,RMED 

UM2,UM3,0M4 

I, CLASS  (I)  ,CFREQ  (I)  ,NGROSS (I) 


C 

C******************************* ********************** 

SUBROUTINE  TKMEAN  (RME AN, R ME A N 1 , RME AN2 , NS) 

COMMON  N , EOH 1  (13)  , PO , B MO, DCS , DTS 

COMMON  FC, FY,ES ,BB,H,DC,DD, AS, AS1 1 ,B1 1 ,D1 1 ,S,C,ZZ, NRU 
COMMON  PHI,EO,J,Z, ECC , EY,FCONCC ( 2 0)  , AS C, EO H  (4 0)  ,FCONST (25) 
COMMON  X  (16  000)  ,  EC  (20)  ,B  (20)  ,P  (4  0)  ,BMM  (1000)  ,BM  (40)  ,FCS  (20) 
COMMON  FST,E  (20)  , N B , DS  (2 0)  , ASB  (2 0)  , FS S  (20 )  , SBM  (20) 

DIMENSION  RMEAN (25) 

C  SET  EACH  VARIABLE  EQUAL  TO  ITS  MEAN  VALUE 
FC=R  ME AN  1 

FY= (RMEAN 2- 4 000. 0) *0.97 
BB  =  RME AN  (3) 

H  =  RMEAN  (4) 

B11  =  RMEAN  (5) 

Dll  =  RMEAN  (6) 

DC=RMEAN  (7) 

DD  =  RMEAN  (8) 

DO  2  1=1, NB 
NI=I+8 

2  DS  (I)  =  R M E A N  (NI) 

C  CALCULATE  THEORETICAL  INTERACTION  DIAGRAM 
CALL  THEORY 

C  WRITE  MEAN  THEORY  INTERACTION  DIAGRAM 
WRITE  (6,100) 
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237 

238 

239 

240 

24  1 

242 

243 

244 

24  5 

246 

247 

246 

249 

250 

251 

252 

253 

254 

255 

256 

2  57 

258 

259 

2oG 

261 

262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 

273 

274 

275 

276 

277 

278 

279 

280 

281 

262 

283 

284 

285 

286 

287 

288 

289 

290 

291 

292 

293 

294 

295 


19C 


100  FORM AT  ('  1  '  ,//////19X,  •  ****ME\N  THEORY  INTERACTION  DIAGRAM****') 
WRITE  (6,514)  NS , NRU 

519  FORMAT  (31X, •  (• , 14,  'SIM' ,15, • ) •//) 

WRITE  (6,101) 

101  FORMAT  (//I  9  X,  '  P  (J)  LBS  '  ,  6  X ,  •  M  (J )  LB-IN' , 7X, • EOH  (J)  ' ) 

DO  4  J  =  1 , N 

WRITE  (6,102)  P  (J)  ,BM  (J)  ,EOH  (J) 

102  FORMAT  (/ 1 6  X  ,  3E  1  5 . 7) 

4  CONTINUE 

WRITE  (6,103) 

103  FORMAT  (//20X,'PO  LBS',7X,'BMO  LB-IN') 

WRITE  (6,104)  PO  ,  B  MO 

104  FORMAT  ( 1 6X , 2E 1 5 . 7) 

RETURN 

END 

C******************************* ********************** 

c 

c***************************************************** 

SUBROUTINE  ACI  (NS) 

C 

C  THIS  SUBROUTINE  CALCULATES  THE  ACI  INTERACTION  DIAGRAM 

C 

COMMON  N,EOH1 (13) , PO , B MO , DCS , DTS 

COMMON  FC,EY,ES ,B3,H, DC,DD,AS, AS  1 1 , B1 1 , D 1 1 , S , C, ZZ , NRU 
COMMON  PHI,EO,J,Z, ECC , E Y , FCONCC  ( 20)  , ASC, F OH (4 0)  , FCONST  (25) 

COMMON  X  (16000)  , EC  (20)  , E  (20) , P  (4 0) , BMM ( 1 000 ) , BM (40 )  ,FCS  (20) 

COMMON  FST , E  (2 0)  , NB , DS  (20 )  , AS B (2 0)  , FSS (2C)  ,SBM  (20) 

IF  (FC.LE.400C .0)  GO  TO  1 
B 1=0. 8 5-0. 05* (FC-4CGC. 0)  /  I  000.0 
IF  (B1.LE.0.65)  B  1  =  0 . 6  5 
GO  TO  4 

1  B 1=0. 8 5 
4  E4=0 .003 

C  CALCULATE  PURE  AXIAL  LOAD  CAPACITY 
PO=0.85*FC*(BB*H-AS) +AS*FY 

C  CALCULATE  AXIAL  LOAD  CAPACITY  AT  BALANCED  CONDITIONS 
PB=0.85*E1*FC*EB*DTS*  (0 . 00  3/ (FY/ES  +  0 . 003) ) 

C  CALCULATE  PURE  MOMENT  CAPACITY 
AST=AS- ASC 
AA=BB*FC*B1*0. 85 
AB=0.003*ASC*ES“AST*FY 
AC=-0. 003*ASC*ES*DCS 
RA=SQPT(AB**2-4 .0*AA*AC)/(2.0*AA) 

C1=(-AB/(2.0*AA))-RA 

IF  (Cl  .  I.E  .0.0)  C  1  =  R  A-  (AB/  (2.  0*AA)  ) 

ES2  =  0. 003*  (Cl- DCS) /Cl 
IF  (ES2.GE. (FY/ES) )  GO  TO  9 

BMO=  ASC *ES 2  *F,S*  (DTS-DCS)  +  (  AST *FY- ASC*ES2  *ES)  *  (DTS- B 1  *C  1/2 . 0) 

GO  TO  8 

9  BMC=  ASC* FY* (DTS-DCS)  +  (  (AST-ASC) *FY) * (DTS- (AST- ASC)  *F Y/ ( 1 . 7 *FC* BE) ) 
C  INITIALIZE  STRAIN  IN  TENSION  STEEL 
8  E 1 =0 . 0 0 19 
J  =0 

2  J=J+1 

IF  (J. EQ. 1)  GO  TO  5 
C  MODIFY  TENSION  STEEL  STRAIN 

IF  (El. GT. -0.001)  E1=E 1-0.0005 
IF  (El  .LF. -0.001)  E1=  E 1- 0 . 00 1 

C  CALCULATE  NEUTRAL  AXIS  DEPTH 
C  =  F4*DD/  (E4-E1) 
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PHI=E9/C 

C  CALCULATE  FORCES  IN  STEEL  BARS 
CAI L  ASTEEL  (E4) 

C  CALCULATE  CONCFETE  COMPRESSIVE  BLOCK  FORCE 
IF  (C.GE. (H/B1) )  C=H/B1 
FCCONC=0.85*FC*B1*BB*C 
C  CALCUI.ATE  AXIAL  LOAD  LEVEL 
P (0) =FCCONC+FST 

C  CALCULATE  BENDING  MOMENT  DUE  TO  CONCFETE  COMPRESSIVE  FORCE 
COMPM=FCCONC*(DD-B1*C/2.0) 

SM  =  0 . 0 

C  CALCULATE  BENDING  MOMENT  DUE  TO  STEEL  FORCES 
DO  3  1=1, NP 

SBM  (I)  =  (FSS  (I)  -FCS  (I)  *  ASB  (I)  )  *  (DD-DS  (I)  ) 

3  SM=SM+SBM (I) 

C  CALCUI.ATE  TOTAL  BENDING  MOMENT  CAPACITY 
BK (0) =C0MPM*SM-P (J) * (DD-H/2. 0) 

GO  TO  6 

5  P (J) =  0 .85*FC*  (  (BB*H)-AS)  ♦  AS*FY 
BM  (J) =0.0 

C  CALCULATE  ECCENTRICITY  E/H 

6  EOK  (J)  =BM  (J)  /  (P  (J)  *H) 

IF  (J.GE.20)  GO  TO  7 

IF  (EOH  (J)  .LT. 2. 0)  GO  TO  2 

7  N= J 

C  WRITE  THE  ACI  INTERACTION  DIAGRAM 
WRITE  (6,100) 

100  FORMAT  (' 1 ' ,//////23X, ' ****ACI  INTERACTION  DIAGRAM****') 
WRITE  (6,519)  NS , NRU 

519  FORMAT  (3  1 X,  •  ( '  ,  14  ,  *  SIM  •  ,15,  • )  •  //) 

WRITE  (6,101) 

101  FORMAT  (//19X,  '  P  (0)  LBS  '  ,  7  X,  '  P  ( B)  LBS  •  ,  6  X  ,  •  M  (0)  LB-IN') 
WRITE  (6,102)  PO , PB, 3  MO 

102  FORMAT  (/I 6X , 3E 1 5 . 7) 

WRITE  (6,103) 

103  FORMAT  (//19X,'P(J)  LBS  '  ,6X, '  M'  (J)  LB-IN  •  ,7X,  •  EOH  (J)  • ) 

DO  20  J= 1 ,  N 

WRITE  (6,104)  P  (J)  ,BM  (J)  ,EOH  (J) 

104  FORMAT  (/ 1 6  X  ,  3E  1  5 . 7) 

20  CONTINUE 

RETURN 

END 

£*********$******************************************* 

SUBROUTINE  ASTEEL  (E4) 

THIS  SUBROUTINE  CALCULATES  THE  ACI  FORCES  IN  THE  STEEL 
COMMON  N,EOH1  (13)  ,  PO , B MO , DCS , DTS 

COMMON  FC , F Y , ES , BB , H , DC , DD , AS , AS  1 1 , B 1 1 , D 1 1 , S , C , ZZ , NF U 
COMMON  PHI,EO,J,Z, ECC , E Y , FCO NCC (20)  , A SC, EOH  (4  0)  ,FCONST  (2  5) 
COMMON  X  (16000)  ,  EC  (2  0)  ,  B  (20)  ,  P  (4 0)  ,  BM M  ( 1  0 00 )  ,  BM  (9 0 )  ,FCS  (20) 
COMMON  FST,E  (20)  ,NB,DS  (20)  , ASB  (2  0) ,FSS  (2  0)  ,SBM  (2  0) 

EY=FY/ES 

FST  =  0 . 0 

DO  4  1=1, NP 

E  (I)  =E4-PHI* DS  (I) 

IF  (DS  (I)  .GE.C)  GO  TO  5 
FCS (I) =0. 85*FC 
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FSC=-FCS  (I)  *  AS  B  (I) 

GO  TO  8 
5  FSC=C . 0 
FCS  (I) =0.0 

8  IF  (E  (I)  .GE.EY)  GO  TO  1 
IF  (E  (I)  .  LE.-EY)  GO  TO  2 
GO  TO  3 

1  FSS  (I)  =FY* ASB  (I) 

GO  TO  4 

2  FSS  (I)  =-FY*  ASB  (I) 

GO  TO  4 

3  FSS  (I)  =E  (I)  *ES*ASB  (I) 

4  FST=FST+FSS  (I)  +FSC 
RETURN 

END 

C********  ********************************************* 

c 

C**  *************************************:t<***  ********** 

SUBROUTINE  PROP  (NS) 

THIS  SUBROUTINE  READS  AND  WRITES  THE  COLUMN  PROPERTIES 
COMMON  N.EOH1  (13)  , PO , B MO , DCS , DTS 

COMMON  FC,FY,ES,BB,H,DC,DD,AS,  AS  1 1 ,  B 1  1 ,  D 1 1  ,  S  ,  C  ,  ZZ  ,  NF.  U 
COMMON  PHI,EC,J,Z,ECC,EY,FCCNCC  (20)  ,ASC,EOH (4  0)  ,FCONST (25) 
COMMON  X  ( 1  6000)  ,EC  (20)  ,  B  (2  0)  ,  P  (4 0)  ,BMM  (1000)  , BM  (40 )  ,FCS  (20) 
COMMON  FST,E  (20)  , NB  ,  DS (20 )  , AS B (2 0)  ,  FSS (2 0)  ,  SBM (20) 

C  FC=CONCRETE  STRENGTH 

C  FY=STEEL  STRENGTH  (PSI) 

C  ES=STEEL  MODULUS  OF  ELASTICITY  (PSI) 

C  BB=CROSS  SECTION  WIDTH  (IN) 

C  H=CROSS  SECTION  DEPTH  (IN) 

C  DC=DISTANCE  TO  THE  COMPRESSION  STEEL  (IN) 

C  DD=DISTANCE  TO  THE  TENSION  STEEL  (IN) 

C  ASC=AREA  OF  COMPRESSION  STEEL  (SQ  IN) 

C  AST= AREA  OF  TENSION  STEEL  (SQ  IN) 

C  AS=TOT AL  AREA  OF  STEEL  (SQ  IN) 

C  AS1 1 =AREA  OF  STIRRUP  (SQ  IN) 

C  B 1 1  =VIIDTH  OF  CORE  (IN) 

C  D11=DEPTH  OF  CORE  (IN) 

C  S=SPACING  OF  STIRRUPS  (IN) 

READ  (5,100)  BB,H,DD,DC,AS,ASC 

100  FORMAT  (6F5.2) 

READ  (5,101)  FC , FY , ES 

101  FORMAT  (3F10.0) 

READ  (5,102)  DCS,DTS,S,B11,D11,AS11 

102  FORMAT  (6F5.2) 

READ  (5,110)  NB 

110  FORMAT  (113) 

DO  1  1=1, NB 

READ  (5,111)  ASB  (I)  ,  DS  (I) 

111  FORMAT  (2F5.2) 

1  CONTINUE 

WRITE  (6,103) 

103  FORMAT  (' 1 1 ,//////3CX ,' COLUMN  CROSS  SECTION  PROPERTIES’) 
WRITE  (6,519)  NS , NRU 

519  FORMAT  ( 3 1 X , ' ( ' , 14 , ' S I M ’ , I 5 , ' ) ' //) 

WRITE  (6,  104)  • 

104  FORMAT  ( 16  X , ’ FC  (PSI)  •  , 3X , • FY  (P SI )  ' , 5 X , • ES  (PSI)’) 

WRITE  (6,105)  FC , F Y , ES 
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105  FOFKAT  (/ 1 6X , 1F6 . 1 , 5X , 1 F7 . 1 , 4 X , 1 F 1 0 . 1) 

WRITE  (6,106) 

106  FORMAT  (//1 6  X  ,  •  B  (IN)  •  ,  5X  ,  •  H  (I  N)  •  ,  SX  ,  '  D  (IN)  •  ,  4X  ,  •  DC  (I  N)  • 

1  , 4X , • AS  (SCIN)  ' , 3X  ,  * ASC  (SQIN)  ' ) 

WRITE  (6,107)  RB,H,DD,DC,AS, ASC 

107  FORMAT  (/ 1  2  X ,  6F 1  0 . 2) 

WRITE  (6,108) 

108  FORMAT  (// 1 6X , *  DCS (IN)  ' , 2 X , ' DT3 ( IN)  ' , 4X, 'S  (IK)  * , 4X , ' B  1  1  (IN)  • , 
14X, • Dll  (IN)  • ,2X,' AS  11  (SQIN) • ) 

WRITE  (6,109)  PCS, DTS, S,B1 1 ,D1 1 , AS1 1 

109  FORMAT  (/I 2X , 6F 1 0 . 2) 

WRITE  (6,112) 

112  FORMAT  (//16X,  •  NB«  ,4X,  •  ASB  (I)  •  ,5X, +  DS  (I)  •) 

DO  1 14  1=1, NF 

WRITE  (6,113)  KB  ,  ASB  (I)  ,  DS  (I) 

113  FORMAT  ( 15X,  113, 2F10. 2) 

114  CONTINUE 
RETURN 
END 


SUBROUTINE  CURVE 


THIS  SUBROUTINE  FITS  A  POLYNOMIAL  TO  THE  INTERACTION 
DIAGRAM 

COMMON  N , EOH 1  (13)  , PO , B MO, DCS , DTS 

COMMON  FC,FY,ES,BB,H,DC,PD,AS, AS11,B11,D1 1 ,S,C,ZZ,NRU 
COMMON  PHI,EO,J,Z,ECC,EY,FCONCC(2  0)  ,ASC,EOH(40)  ,FCONST  (25) 
COMMON  X  (16000)  ,EC(20)  ,B(20)  ,P(40)  ,BMM  (1  000)  ,  BM  (40)  ,FCS  (20) 
COMMON  FST , E  (20)  ,NB,DS  (20)  ,ASB  (20)  ,FSS  (20)  ,SEM  (20) 

DIMENSION  DI  (4  00)  ,D(7C)  ,  SB  (10)  ,T  (10)  ,COE  (11) 

DIMENSION  XBAR  (11)  , STD  (11)  , SUMSQ  ( 1 1 ) , IS A VE ( 1 1 )  , A NS  (10) 
DIMENSION  XX  (500)  ,BBB  (10) , EE  (10) 

M  =  1  0 

EOH 1  (1) =0.05 
EOH 1 (2) =0.  10 
EOH 1  (3) =0.15 
EOH 1  (4) =0. 20 
EOH 1  (5) =0.30 
EOH1  (6) =0. 40 
F.OK1  (7)  =0.50 
EOH 1  (8) =0.60 
EOH 1  (9) =0.70 
EOH  1  (10) =0.80 
EOH 1  (11) =0 .90 
EOH 1  (12) =1.00 
EOH  1  (13)  =1 . 50 
N  =  N-  1 

DO  2  J  =  1  ,  N 
P(J)  =P  (J+1) 

BM  (J)  =BM  (J+1) 

2  EOH  (J)  =EOH  (J+1) 

DO  10C  1=1, N 
D I F il  =  B M  (I  +  1)  -BM  (I) 

IF(DIFM)  105,100,100 
100  CONTINUE 
105  N  P=I 

EOHE  =  ECH  (I) 


476 

477 

478 

479 

480 

481 

482 

483 

484 

485 

486 

467 

488 

489 

490 

491 

492 

493 

494 

495 

496 

497 

498 

499 

500 

501 

502 

503 

504 

505 

506 

507 

508 

509 

510 

511 

512 

513 

514 

515 

516 

517 

518 

519 

520 

521 

522 

523 

524 

525 

526 

527 

528 

529 

530 

531 

532 

533 

534 

535 


194 


mi  =  m 


IF  (M. GE.  (NP- 1) )  M  =  NP-2 
L=NP*n 

DO  110  1=1, NP 
J  =  L  ♦  I 

XX  (I)  =  EOH  (I) 

110  XX(J)=P(I) 

LL=L+NP 

CALL  GDATA  ( NP , H , XX , XB AR , STD , D , SUMSQ) 

MM  =  M  -*■  1 


SUM  =  0. 0 
DO  200  1=1, n 
ISAVE  (I)  =1 

CALL  ORDER  (  HP!,  D  ,  MM,  I  ,  IS  AV  E,  DI,  EE) 

CALL  H-INV  (DI,  I,  DET,BBB,T) 

CALL  TMULTR  ( N P , I , XB AR , STD , SUM SQ , DI , EE , IS A VE , BBB , S B , T , A NS ) 
IF  (ANS  (7) )  220,130,130 
130  SUHIP  =  ANS  (4) -SUM 

IF  (SOMIP)  220,220,150 
150  SUM  =  ANS  (4) 

COE  (1)  =ANS  (1) 

DO  160  J  =  1,I 
160  COE  (J+1)  =BBB  (J) 

11=14-1 

JJ=I41 

200  CONTINUE 
220  NN= 1 3 

DO  240  11=1, NN 
ECH2=EOHl  (II) 

IF  (EOH2.GT.EOHB)  GO  TO  250 
P  (II) =0.0 
DO  245  1=1, JJ 

245  P  (II)  =P  (II)  +CO E  (I)  *  (EOH1  (II)  **  (1-1)  ) 

240  CONTINUE 
250  NNN=II 

DO  101  I=NP , N 
DIFE=3 . 0-EOH  (I) 

IF  (DIFE)  102,101,101 

101  CONTINUE 
GO  TO  103 

102  N=I— 1 

103  NM=N-NP4l 
M  =  M  1 

IF  (M.GE.  (NM-1) )  M=NM- 2 
L=NH*M 

DO  310  1=1, NM 
J  =  L  +  I 

XX (I) =1 .0/EOH  (I  +  NP-1) 

310  XX (J) =BM(I+NP-1) 

CALL  GDATA  (N M , M , XX , XB AR , STD, D , S UMSQ) 

M  M  =  M ♦  1 
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SUM=0 . 0 
DO  300  1=1, M 
ISAVE  (I)  =1 

CALL  ORDER  (MM , D , MM , I , IS AVE, DI , EE) 
pa  t  T  MTNV  ET / B  E  B / T ) 

CALL  ?SSlTR  (N-M.I,XBAR,STD,SUMSQ,DI,EE,ISAVE,BBB,SB,T,ANS) 

IF  (ANS  (7 ) )  320,330,  330 

SUMI P=  ANS (4) -SUM 
IF  (SUMIP)  320,320,350 
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350  SUM=ANS(4) 

COE  (1) =  ANS  (1) 

DO  360  J  =  1,I 
360  COE  (J+  1)  =BBB  (J) 

JJ=I+1 

300  CONTINUE 
320  DO  340  II=NNN, NN 
BM  (II) =0.0 
DO  345  1=1, JO 

345  BH  (II)  =BM  (II)  +  COE  (I)  *  (  (1 .0/EOH1  (II) ) **  (I- 1) ) 

P  (II)  =  BM  (II)  /  ( H*  EOH 1  (II)) 

340  CONTINUE 
N4=NNN— 1 
DO  370  1=1, N4 
370  BM(I)=P(I)  *H*EOH 1  (I) 

RETURN 

END 

C* ****************************************** ********** 

C 

C**************************»************************** 

SUBROUTINE  RANDOM  (IY , SD , SB, CONST, ITP,  V) 

C 

C  THIS  SUBROUTINE  GENERATES  VALUE  OF  THE  VARIABLES  KITH 
C  THE  MEAN,  STANDARD  DEVIATION,  AND  DISTRIBUTION  GIVEN 

C 

COMMON  N,EOHl  (13)  ,PO,BMO, DCS, DTS 

COMMON  FC,FY,ES ,BB,H, DC,DD,AS, AS  11 ,B1 1 ,D1 1 ,S,C,ZZ,NRU 
COMMON  PHI , EO , J , Z , ECC , EY,FCONCC (20) ,ASC,EOH(40) ,FCONST(25) 
COMMON  X  (16C00)  ,  EC  (20)  ,B  (20)  ,P  (4  0)  ,BMM  (1000)  ,BM  (40)  ,FCS  (20) 
COMMON  FST , E  (20)  ,NB,DS  (20)  ,ASB  (20) ,FSS (20)  ,SBM (20) 

A=0.0 

DO  50  1=1,12 
IY=I Y*65539 

IF  (IY.LT.0)  IY=IY+2147483647+  1 

Y=IY*0.4656613E-9 
50  A= A  ♦  Y 
A = A  -  6 . 0 
V=A*SD+RM 

IF  (ITP. EQ . 1 )  V=  1 0 . 0**  V 

IF  (ITP.EQ.2)  V=1 0. 0**V+CONST 

RETURN 

END 

C* ******************** *********************  *********** 

C 

£***************************************************** 

SUBROUTINE  THEORY 


C 

c 

c 


c 


THIS  SUBROUTINE  CALCULATES  THE  THEORETICAL  P-M  DIAGRAM 


N, EOH  1(13)  ,PO,BMO,DCS,DTS 

FC ,FY,ES,BB,H,DC,DD,AS,AS11,B11,D11,S,C,ZZ,NRU 

PHI , EO , J , Z , ECC ,EY,FCONCC  (20)  ,ASC,EOH(4  0)  ,FCONST  (2  5) 

X (16000)  , EC  (20)  ,B  (20) , P (40) , BMM (1000) , BM (40)  ,FCS  (20) 

FST,  E  (20)  ,  NB,  DS  (2  0)  ,  AS B  (2 0)  ,  FSS  (2 0 )  ,SBP1  (20) 

FC*0. 85 

(FC.EQ. 1000.0)  FC= 1000. 1 


COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
FC 
IF 

ECC=5700C.0*SQRT  (FC) 


EO= 1 . 8  *FC/ECC 
EY=FY/ES 

CALCULATE  PUPE  AXIAL  LOAD  CAPACITY 
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PO=FC*BE*H+AS*FY-AS*FC 

J=0 

1  J=J*1 

C  SET  AXIAL  LOAD  LEVEL 

IF  (J.  EQ.  1)  GO  TO  5 
IF  (J.EQ.21)  P  ( J)  =0 . 0 

IF  (P  (J-1)  .LE.  (0.6*PO)  )  P(J)=P(J-1)-0.034*PO 
IF  (  P  ( J- 1 )  .LE.  (C.  1  *PO)  )  P(J)  =  P  (J  - 1)  -0  .  04  *  PC 
IF  (P  (J-  1)  ,GT.  (0.6  +  PO)  )  P  (J)  =P  (J-  1) -0.  16*PO 
IF  (P  (J)  .LE.0.  C)  P  ( J)  =0.0 
IF  (J.EQ.2)  P(J)  =P(J-1)-0.08*PO 
C  CALCULATE  MOMENT  CAPACITY  AT  SPECIFIED  AXIAL  LOAD 
CALL  AXIAL 

IF  (P  (J)  . EQ. 0. 0)  BMO=BM  (J) 

IF  (P(J)  .  EQ .  0 . 0)  GO  TO  7 
GO  TO  6 

5  P(J)=PO 
BM  (J) =0.0 

6  EOH  (O)  =BM  (J)/(P  (J)  *H) 

GO  TO  1 

7  N=J-1 

C  ELIMINATE  EREATIC  POINTS  ON  THE  INTERACTION  CURVE 
M  =  N 

HJ  =  N-  1 

DO  8  I  J=3 , N J 

IF  (BM  (IJ)  .GE.BM  (IJ-1)  )  GO  TO  8 
IF  (BM  (IJ+  1)  .LE.BM  (IJ-  1)  )  GO  TO  8 
M  =  M-1 

DO  9  J JJ=I J ,  N.I 
P (JJJ) =P  (JJJ  +  1) 

BM  (JJJ)  =BM  (JJJ+  1) 

9  EOH  (JJJ)  =EOH  (JJJ*  1) 

8  CONTINUE 
N  =  M 

RETURN 

END 

C******** $  A******************************* ********** 

***************************************************** 

SUBROUTINE  AXIAL 

THIS  SUBROUTINE  CALCULATES  THE  MOMENT  AFTER  BALANCING  P 
COMMON  N,EOH1 (13)  , PO , B MO , DCS , DTS 

COMMON  FC,FY,ES ,BB,H, DC, DD, AS, AS  11,31 1,D1 1 ,S,C,ZZ,  NRU 
COMMON  PHI,EO, J,Z, ECC, EY.FCONCC (20) ,ASC,EOH(40) ,FCONST(25) 
COMMON  X  (16000)  , EC  (2  0)  ,B  (2  0)  , P  (4  0)  , 3MM (1000)  , BM  (4  0)  , FCS  (20) 
COMMON  FST,E  (2  0)  ,NB,DS  (20)  ,ASB  (2  0)  ,FSS  (20)  ,SBM  (20) 

PK 1=0.000000 1 
PHIH=PHI*H 
11=  1 

14  E4=0 .002 
EINCR=0.002 
33  E4  =  E4- El  NCR 

El NCR=EIM CP/2.0 
32  E4  =  E4  +  EINCF( 

,FCCOffC  =  0. 0 
C=F4 /PHI 
ECO=  (C-H) *  PHI 
IF  (C.GE.H)  C=H 
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IF  (C.LT.H)  ECO=0 . 0 
ASC=0. 0 
DO  34  1=1, NB 

IF  (DS(I).LE.C)  ASC=ASC*ASB (I) 

34  CONTINUE 

C  CALCULATE  PARAMETERS  OF  THE  CONCRETE  STRESS  STRAIN  CURVE 
P11=(2.0*(B11«-D11)  *  AS  1  1  ♦  ASC*S)  /(B11*D11*S) 
E5CH=0.75*P11*SQRT (B11/S) 

E50U=(3. 0+0. 002*FC)/(FC- 1000.0) 

IF  (E50U.LE.0.0)  E50U=0.06 

Z=0.5/(E50H«-E50U-EO) 

ZZ=0.5/(E50U-EO) 

EU  =  2. 0*  (E50H+E50U) -EO 
DX=C/10 

C  CALCULATE  THE  CONCRETE  COMPRESSION  BLOCK  FORCE 
DO  23  1=1,10 
AI  =  I 

X  (I)  =C- AI *  DX  +  DX/2 
EC  (I)  =  PHI*  X  (I)  «•  ECO 
B  (I)  =BB 

C  MAXIMUM  STRAIN  FOR  UNCONFINED  COMPRESSION  0.004 
IF  (EC  (I)  .LE.EO)  GO  TO  3 
IF  (EC  (I)  .GE. 0.004)  GO  TO  21 
IF  (EC  (I)  .GT.EO)  GO  TO  4 

3  FCC=  FC * (2. 0*EC (I) /EO- (EC  (I) /EO) **2) 

GO  TO  22 

4  FCC=FC*  (1 .  0-  Z*  (EC  (I)-EO)  ) 

ECU  =  FC*  (1 . 0-ZZ*  (EC  (I) -EO) ) 

IF  (FCC.LE.0.0)  FCC=0 . 0 
IF  (FCU.LE.0.0)  FCU=0 . 0 

FCONCC  (I)  =FCC*DX*B11*FCU*DX*  (B  (I)  -B  1 1) 

GO  TO  23 

21  B  (I)  =B  1 1 

FCC=FC*  (1 . 0-Z*  (EC  (I)  -EO)  ) 

IF  (FCC.LE.0.0)  FCC=0 . 0 
IF  (X  (I)  .GE.  (C-DC)  )  ?CC=0.0 

22  FCONCC  (I) =FCC*DX*B (I) 

23  FCCONC=FCCCNC+FCONCC (I) 

C  CALCULATE  THE  CONCRETE  TENSION  BLOCK  FORCE 
IF  (C.GE.H)  GO  TO  25 
SFC  =  SQRT  (FC) 

ET=7 . 5*SFC/ECC 
TC=ET/PHI 
TCA=H-C 
RTC  =  TC  A/TC 

IF  (TC.GT.TCA)  TC=TCA 
IF  (TC.LE.TCA)  RTC=1.0 
FCONCT=-RTC*7 . 5*SFC*TC/2 . 0*BB 
GO  TO  18 
25  FCONCT=0 . 0 
18  CALL  FSTEEL  (E4) 

C  CHECK  FORCE  COMPATIBILITY 

PAXIAL=FCCONC*FCONCT+FST 

TOLA  =  P  (J) *0.02 

IF  (P (J)  .EQ.0.0)  TOLA=  0.00 1* PO 
TOL=P  (.7)  -PAXI  AL 
IF  (TOL.LT.-TOLA)  GO  TO  33 
IF  (TOL.GT.TOLA)  GO  TO  35 
GO  TO  36 

35  IF  (F4.GE.EU)  GO  TO  44 
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IF (EINCR . GE. 0  . 0C00001)  GO  TO  32 

36  compm=o.o 

C  CALCULATE  THE  MOMENT  DUE  TO  THE  CONCRETE  COMPRESSION  FORCE 
DO  24  1=1,10 

24  COMP  M  =  COM PM ♦ FCONCC (I) *  (DD-C+X (I) ) 

SM  =  0 . 0 

C  CALCULATE  THE  MOMENT  DUE  TO  THE  STEEL  FORCES 
DO  13  I=1,NB 

SBM (I)  =  (FSS (I) - FCS  (I)  ♦ A  SB  (I) )  * (DD-DS (I) ) 

13  SM=SM+  SBM  (I) 

IF  (FCONCT.EQ.0.0)  TC=0.0 
C  SUM  THE  MOMENTS  ABOUT  THE  TENSION  STEEL 

BMM (II) =COMPM+SM+FCONCT*  (DD-  (C+2.0+TC/3. 0) )-P  (J) * (DD-H/2.0) 
IF  (II. EQ. 1)  GO  TO  17 
TOLBMA  =  ABS  (BMM  (II-  1) *0 . 01) 

BMTOL=BMM  (II)-BMM  (II-1) 

IF  (BMM(II) .LE.0.0)  GO  TO  42 
IF  (BMTOL.GE.  (0. 5*TOLBMA) )  GO  TO  41 
IF  (BMTOL. LE.-TOLBMA)  GO  TO  42 
GO  TO  16 
17  PHI NCR=0 .001 
GO  TO  41 
44  E4=0 .001 

PHIH=PHIH— PHINCR 
PHI NCR=PHI NCR/5 . 0 
PHIH=PHIH+PHINCR 
PHI=PHIH/H 
El NCR=EI NCR/2 . 0 
GO  TO  32 

42  PHI H= PHI H- PHINCR 
PHI  NCR  =  PHI NCR/5. 0 
41  PHIH=Ph'IH  +  PHINCR 
PHI=PHIH/H 
11=11+1 
GO  TO  14 

16  BM  (J)  =BMM  (II- 1) 

CONTINUE 

RETURN 

END 

Q&  &$#*#  + $*******#********  ***********************  £#*■$*# 

c 

SUBROUTINE  FSTEEL  (E4) 

C 

C  THIS  SUBROUTINE  CALCULATES  THE  THEORY  FORCES  IN  THE  STEEL 

C 

COMMON  N,EOHl  (13)  , PO ,BMO,DCS,DTS 

COMMON  FC,FY,ES , BB, H, DC, DD, AS, AS1 1 ,B1 1 ,D1 1 ,S,C,ZZ, NRU 
COMMON  PHI, EO, J, Z, ECC, EY, FCONCC (20) ,ASC,EOH(40) ,FCONST(25) 
COMMON  X (16000)  , EC  (20)  , B  (2 0)  , P  (4 0) , BM M  ( 1 00 0)  , BM  (4  0 )  , FCS  (20) 
COMMON  FST,E  (20)  ,NB,DS  (20)  ,ASB  (2  0)  ,FSS  (2  0)  ,SEM  (20) 

FST=0 . 0 

DO  4  1=1, NB 

E (I) =E4-PHI*DS (I) 

IF  (DS  (I)  .GE.C)  GO  TO  5 
IF  (E  (I)  .GT.  EO)  GO  TO  6 
FCS (I) =  FC*  (2 . 0*E  (I) /EO- (E (I)/EO)  **2) 

GO  TO  7 

6  FCS (I) =FC* (1.0-Z* (E  (I) -EO) ) 

IF  (FCS  (I)  .  LE.  0. 0)  FCS  (I)  =0.0 
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7  FSC=-FCS  (I)  *ASB  (I) 

GO  TO  3 

5  FSC=0 . 0 
FCS (I) =0.0 

8  E  (I )  =-  F,  (I) 

IF  (E (I) . GE. EY)  GO  TO  2 
IF  (E  (I)  .  LF.-EY)  GO  TO  1 
GO  TO  3 

1  FSS  (I) =FY*ASB (I) 

GO  TO  4 

2  FSS  (I) =-FY*ASB (I) 

GO  TO  4 

3  FSS  (I)  =-E  (I)  *ES* ASB  (I) 

4  FST=FST  +  FSS  (I)  +FSC 
RETURN 

END 

C+*******:*‘*****’**,$'*<'********;<c*************’$l’tt>!'<‘*****;fc,l‘I!< 

c 

C**********;*i**'S'*******’<‘********************;<'******#*** 

SUBROUTINE  ST AT  (F a , ST D V ,COV , COS ,COK, R MI N , EM AX , UM2 , UM3 , U H4 , IM AX , I M 
1IN, BMED,CLASS, CFFEQ, NG  ROSS) 

C 

C  THIS  SUBROUTINE  CALCULATES  THE  MEAN,  COEFFICIENT  OF 
C  VARIATION,  COEFFICIENT  CF  SKENNESS,  COEFFICIENT  OF 
C  KUPTOSTS, AND  CUMULATIVE  FREQUENCY  TABLE 
C 

COMMON  N , EOH  1  (13)  , PO , BMO , DCS , DTS 

COMMON  FC,FY,ES,BB,H,DC,DD, AS, AS11,B11,D11 ,S,C,ZZ, NFU 
COMMON  PHI,EO,J,Z,ECC, EY ,FCONCC (20)  ,ASC,EOH(40)  ,FCONST  (25) 

COMMON  X  (16000)  , EC  (20)  ,B  (2  0)  , P  (4 C) , EMM ( 1 000) , BM (40)  ,FCS  (20) 

COMMON  FST,E  (2  0)  ,NB,DS  (20)  ,ASB  (2  0)  ,FSS  (2  0)  ,SBM  (20) 

DIMENSION  CLASS  (31)  ,CFREQ  (31) 

DIMENSION  Y  (800  2)  , DIF (1600  0) , NGFOSS (31) 

SUM=0 . 0 
DO  10  1=1, N 
10  SUM=  SUM  +  X (I) 

RM=SUM/N 
S  UM  =  0 . 0 
DO  20  1=1, N 

20  SUK=SUM+ (X (I)-RM) **2 
U  M2  =  SUM/N 

STDV  =  SQPT  (SUM/  (N-  1)  ) 

CO V=  STD V/R  M 
SUM=C . 0 
DO  30  1=1, N 

30  SUM=SUM+  (X  (I)-RM)  **3 
UM3= SUM/N 

COS  =  SUM/(N*  (ST  DV  **  3)  ) 

S0M=0. 0 
DO  40  1=1, N 

40  SUM=SUM+  (X  (I)-RM)  **4 
UM4  =  SUM/N 

C OK  =  SU M/  (N *  (STDVv  +  U) ) 

RM AX  =  X  (1) 

IMA  X= 1 
DO  50  1=2, N 

IF  (RMAX.  T.T.  X  (I)  )  IMAX  =  I 
50  RM AX  =  X  (I MAX) 

R MI N  =  X  (1) 

IMI N  = 1 


836 

837 

838 

839 

840 

841 

842 

843 

844 

845 

846 

847 

848 

849 

850 

851 

852 

853 

854 

8  55 

856 

857 

858 

8  59 

860 

861 

862 

863 

864 

865 

866 

867 

868 

869 

870 

871 

87  2 

873 

874 

875 

877 

878 

879 

880 

881 

882 

883 

884 

885 

886 

887 

888 

889 

890 

891 

892 

893 

894 

895 

896 


200 


DO  60  1  =  2, N 

IF  (RMIN.GT.  X  (I)  )  IMIN=I 
60  R  M I N  =  X (IM I N ) 

Cl NT=0 . 05 
CLASS  (1) =0.50 
DO  70  1=2,30 

70  CLASS  (I)  =CLASS  (I-  1)  ♦CINT 
CLASS (31) =2.0 
DO  90  11=1,31 
NPOI NT  =  0 
DO  80  1=1, N 

80  IF  (CLASS  (II)  .GT.  X  (I)  )  NPOINT=  NPOINT-*- 1 
NGROSS  (II) = NPOI NT 
90  CFEEQ  (II)  =  (100. 0*NPOINT) /N 
DO  100  1=1,31 

100  CLASS  (I) =CLASS (I) -0.00001 

Y  (1)  =PBIN 
N 1=N/2+2 

DO  66  11=2, N1 
NPOI NT=1 

DIF  (1)  =X  ( 1 )  - Y  (II-  1) 

DO  65  1=2, N 

IF  (X  (I)  . EQ. 0. 0)  GO  TO  65 
DIF  (I)  =X  (I)-Y(II-I) 

IF  (DIF  (I)  .LT. DIF (NPOINT) )  NPOINT=I 

65  CONTINUE 

Y  (II) =X  (NPOINT) 

X (NPOINT) =0.0 

66  IF  (NPOINT. EQ. 1)  X  (NPOINT) =R MAX* 1 . 0 
N  1  =  N  1-1 

DO  67  1= 1 , N 1 

67  Y  (I )  =  Y  (I*  1 ) 

RMED  =  (Y  ( N 1 —  1 ) +Y ( N 1 ) ) /2 . 0 
RN  =  (N/2. 0) - (N/2) 

IF  (RN .  GT .0.1)  RMED=Y(N1) 

RETURN 

END 

C  SUBROUTINE  MULTR 

C*  **********************************************  ****** 

c 

C  PURPOSE 

C  PERFORM  A  MULTIPLE  LINEAR  REGRESSION  ANALYSIS  FOR  A 

C  DEPENDENT  VARIABLE  AND  A  SET  OF  INDEPENDENT  VARIABLES.  THIS 

C  SUBROUTINE  IS  NORMALLY  USED  IN  THE  PERFORMANCE  CF  MULTIPLE 

C  AND  POLYNOMIAL  REGRESSION  ANALYSES. 

C 

C  USAGE 

C  CALL  MULTR  (N , K , X B AR ,STD , D , RX , RY , I S A VE, B , SB , T , A NS ) 

C 

C  DESCRIPTION  OF  PARAMETERS 

C  N 

C  K 

C  XBAR 

C 

C  STD 

C 

C  D 

C 
C 


-  NUMBER  OF  OBSERVATIONS. 

-  NUMBER  OF  INDEPENDENT  VARIABLES  IN  THIS  REGRESSION. 

-  INPUT  VECTOR  OF  LENGTH  M  CONTAINING  MEANS  OF  ALL 

VARIABLES.  H  IS  NUMBER  OF  VARIABLES  IN  OBSERVATIONS. 

-  INPUT  VECTOR  OF  LENGTH  M  CONTAINING  STANDARD  DEVI¬ 

ATIONS  OF  ALL  VARIABLES. 

-  INPUT  VECTOR  OF  LENGTH  M  CONTAINING  THE  DIAGONAL  OF 
THE  MATRIX  OF  SUMS  OF  CROSS-PRODUCTS  OF  DEVIATIONS 
FROM  MEANS  FOR  ALL  VARIABLES. 
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c 

c 

c 

c 

c 

c 

c 

c 


RX  -  INPUT  MATRIX  (K  X  K)  CONTAINING  THE  I 
INTERCORRELATIONS  AMONG  INDEPENDENT  V 
RT  -  INPUT  VECTOR  OF  LENGTH  K  CONTAINING  I 
TIONS  OF  INDEPENDENT  VARIABLES  WITH  D 
VARIABLE. 

ISAVE  -  INPUT  VECTOR  OF  LENGTH  K+1  CONTAINING 
INDEPENDENT  VARIABLES  IN  ASCENDING  OP 
SUBSCRIPT  OF  THE  DEPENDENT  VARIABLE  I 
THE  LAST,  K  «■  1 ,  POSITION. 

B  -  OUTPUT  VECTOR  OF  LENGTH  K  CONTAINING 

COEFFICIENTS . 

SB  -  OUTPUT  VECTOR  OF  LENGTH  K  CONTAINING 
DEVIATIONS  OF  REGRESSION  COEFFICIENTS 
T  -  OUTPUT  VECTOR  OF  LENGTH  K  CONTAINING 

ANS  -  OUTPUT  VECTOR  OF  LENGTH  10  CONTAINING 

INFORMATION. . 

ANS  ( 1)  INTERCEPT 

ANS  (2)  MULTIPLE  CORRELATION  CCEFFICI 

ANS  (3)  STANDARD  ERROR  OF  ESTIMATE 

ANS  (9)  SUM  OF  SQUARES  ATTRIBUTABLE  T 
SION  (SSAR) 

ANS  (5)  DEGREES  OF  FREEDOM  ASSOCIATED 
ANS  (6)  MEAN  SQUARE  OF  SSAR 
ANS  (7)  SUM  OF  SQUARES  OF  DEVIATIONS 
SION  (SSDR) 

ANS  (8)  DEGREES  OF  FREEDOM  ASSOCIATED 
ANS  (9)  MEAN  SQUARE  OF  SSDR 
ANS  ( 10)  F-VALUE 

REMARKS 

N  MUST  BE  GREATER  THAN  K+ 1 . 


NVFKSF  OF 
API ABLES. 
NTERCORKELA- 
EPENDENT 

SUBSCRIPTS  OF 
DER.  THE 
S  STORED  IN 

REGRESSION 

STANDARD 

T-VALUES. 

THE  FOLLOWING 


ENT 

O  REGRES- 
WITH  SSAR 
FROM  REGRES- 
WITH  SSDR 


SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 
NONE 


METHOD 

THE  GAUSS-JORDAN  METHOD  IS  USED  IN  THE  SOLUTION  OF  THE 
NORMAL  EQUATIONS.  REFER  TO  W.  W.  COOLEY  AND  P.  R.  LOHNES, 
'MULTIVARIATE  PROCEDURES  FOR  THE  BEHAVIORAL  SCIENCES', 

JOHN  WILEY  AND  SONS,  1962,  CHAPTER  3,  AND  B.  OSTLE, 
•STATISTICS  IN  RESEARCH',  THE  IOWA  STATE  COLLEGE  PRESS, 
1959,  CHAPTER  8. 


SUBROUTINE  TMULTR  {N PN , K, X BAR , STD, D , P X , R Y , IS A VE , BB B , S B , T , A  NS ) 
COMMON  N,EOHl (13) ,PO,BMO, DCS , DTS 

COMMON  FC, FY,ES,BB,H,DC,DD,AS,AS11,B11,D11,S,C,ZZ,NFU 
COMMON  PHI , EO , J , Z , ECC, EY , FCONCC ( 20) ,ASC,EOH(40) ,FCONST(25) 
COMMON  X  (16000)  ,EC  (20)  ,B  (20)  ,P  (40)  ,BMM  (1000)  ,BM(U0)  ,FCS(20) 
COMMON  FST ,  E  (20)  ,N3,DS(20)  ,ASB(20)  ,FSS(2C)  ,SBM(2C) 

DIMENSION  XBAE  (11)  ,STD  (1 1)  ,D(11)  ,  RX  (4C0)  ,  R  Y  (1C) 

DIMENSION  ISAVE(II)  ,EBB(10),SB(1C)  ,  T  ( 1  0 )  ,  ANS  ( 10) 


IF  A  DOUBLE  PRECISION  VERSION  OF  THIS  ROUTINE  IS  DESIRED,  THE 
C  IN  COLUMN  1  SHOULD  BE  REMOVED  FROM  THE  DOUBLE  PRECISION 
STATEMENT  WHICH  FOLLOWS. 


DOUBLE  PRECISION  XE A P , ST D , D, R X , R Y , B , SB , T , 
FN,FK,SSARH,SSDRM,F 


ANS, RM, BO, SSAR, SSDR, SY, 


956 

959 

960 

961 

962 

963 

964 

965 

967 

968 

969 

971 

972 

973 

979 

975 

976 

977 

978 

979 

980 

961 

982 

983 

984 

985 

986 

987 

988 

989 

990 

991 

992 

993 

994 

995 

996 

997 

998 

999 

1000 

1001 

1002 

1003 

1004 

1005 

1006 

1007 

1008 

1009 

1010 

1011 

1012 

1013 

1014 

1015 

1016 

1017 

1018 

1019 


20  2 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


1  00 


110 


c 

c 

c 


c 

c 

c 


c 

c 

c 

120 

c 

c 

c 


THE  C  HOST  ALSO  BE  REMOVED  FF  OK  DOUBLE  PRECISION  STATEMENTS 
APPFARING  IN  OTHER  ROUTINES  USED  IN  CONJUNCTION  WITH  THIS 
ROUTINE. 

THE  DOUBLE  PRECISION  VERSION  OF  THIS  SUBROUTINE  MUST  ALSO 
CONTAIN  DOUBLE  PRECISION  FORTRAN  FUNCTIONS.  SORT  AND  ABS  IN 
STATEMENTS  122,  125,  AND  135  MUST  BE  CHANGED  TO  DSQRT  AND  DABS. 


KM  =  K ♦  1 

BETA  WEIGHTS 

DO  100  J=  1  ,  K 
BBB  (J) =0.0 
DO  110  J=  1 , K 
L1  =  K*  (J-  1) 

DO  110  1=1, K 
L=I.  1+1 

BBE  (J)  =B3B  (J)  +RY  (I)  *FX  (L) 

RM  =  0 . 0 
BO=0 .0 
L1=I SAVE  (MM) 

COEFFICIENT  OF  DETERMINATION 

CO  120  1=1, K 

RM  =  RM  +  BBB  (I)  *RY  (I) 

REGRESSION  COEFFICIENTS 

L  =  ISAVE  (I) 

BBB  (I)  =BBB  (I)  *  (STD  (LI)  /STD  (L)  ) 

INTERCEPT 

BO=EO+BBB (I) *XBAR (L) 

BO=XBAR (Li) -BO 

SUM  OF  SQUARES  ATTRIBUTABLE  TO  REGRESSION 
SSAR  =  RM*D  (LI) 


C 

c 

c 

122 


MULTIPLE  CORRELATION  COEFFICIENT 
RM=  SQRT  (  ARS(RK)) 


SUM  OF  SQUARES  OF  DEVIATIONS  FROM  REGRESSION 

SSDR  =  D  (LI)  -SSAR 
IF  (SSDR.EQ.0.0)  SSBR=0.1 

VARIANCE  OF  ESTIMATE 

FN  =  N  PN-K- 1 
SY  =  S  SDR/FN 


C  STANDARD  DEVIATIONS  OF  REGRESSION  COEFFICIENTS 

C 

DO  130  J= 1 »  K 


203 


1020 

L 1  =  K  *  ( J~  1 )  ♦  J 

1  021 

L=I SAVE  (J) 

1022 

125 

SB  (J)  =  SQRT  (  A 

1023 

C 

1024 

c 

COMPUTED  T- 

1025 

c 

10  26 

130 

T  (J)  =BBB  (J)  /SB 

1028 

c 

STANDARD  ER 

1030 

135 

S  Y  =  SQRT (  ABS ( 

1032 

c 

F  VALUE 

1034 

FK=K 

1035 

SSARM=SSAR/FK 

1036 

SSDRM=SSDR/FN 

1037 

F=SS ARM/S SDR M 

1039 

ANS  (1)  =BO 

1040 

ANS  (2) =RM 

1041 

ANS  (3)  =SY 

1042 

ANS  (4) =SSAR 

1043 

ANS  (5) =FK 

1044 

ANS  (6) =SSARM 

1045 

ANS  (7) =SSDR 

1046 

ANS (8) =FN 

1047 

ANS  (9) -SSDRM 

1048 

ANS  (10) =F 

1  049 

RETURN 

1050 

END 

OF  ESTIMATE 


END  OF  FILE 
$SIGNOFF 


appendix  e 

DATA 

INPUT  FOR  THE  MONTE  CARLO  PROGRAM 

Note : 

All  units  are  in  inches  and  pounds. 

Card 

Columns 

Data  Description 

Format 

1 

1-  5 

Number  of  Variables  (NV) 

15 

6-10 

Number  of  Simulations  (NS) 

15 

1 1-19 

Limiting  Steel  Strength  (FY1) 

F9.2 

Note: 

The  limiting  steel  strength  is  a  maximum 

i  value  of 

steel  strength  which  could  reasonably  be  expected. 

This  is 

required  so  that  extremely  high 

values  of 

steel  strength  are  not  used  for  the 

theoretical 

calculations. 

2 

1-15 

Mean  Concrete  Strength  (RMEAN1) 

F15.5 

15-30 

Mean  Steel  Strength  (RMEAN2) 

F15.5 

31-40 

Initial  Seed  (IY) 

110 

41-50 

Number  of  Fun  (NRU) 

110 

Note : 

The  initial  seed  is  any  integer.  This 

number  is 

required 

to  initiate  the  random  number 

generating 

subroutine.  The  number  of  run  is  any 

identifying 

number  for  the  specific  run. 

3 

1-  5 

Width  of  Column  (BB) 

F5 . 2 

6- 1  0 

Depth  of  Column  (H) 

F5. 2 

11-15 

Distance  From  the  Compression  Face 

to 

Longitudinal  Steel  Closest  to  the 

Tension  Face  (DD) 

F5.  2 

16-20 

Distance  From  Compression  Face  to 

Nearest  Longitudinal  Steel  (DC) 


F5.2 
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205 


DATA  INPUT  CONTINUED 


Card 


6 

7 


Note 


8 


Columns 
2 1-25 
26-30 

1-10 
1 1-20 
21-30 
1-  5 

6-10 


Data  Description 

Total  Longitudinal  Steel  Area  (AS) 

Longitudinal  Compression  Steel  Area 
(ASC) 

Concrete  Design  Strength  (FC) 

Steel  Yield  Strength  (FY) 

Steel  Modulus  of  Elasticity  (ES) 

Depth  From  Compression  Face  to  the 
Centroid  of  Compression  Steel  (DCS) 
Depth  From  Compression  Face  to  the 
Centroid  of  Tension  Steel  (DTS) 

11-15  Spacing  of  Steel  Ties  (S) 

16-20  Width  of  Ties  (B11) 

21-25  Depth  of  Ties  (Dll) 

26-30  Area  of  Steel  Tie  (AS11) 

Number  of  Longitudinal  Bars  (NB) 

Area  of  Individual  Steel  Bars  (ASB(I)) 
6-11  Distance  From  Compression  Face  to  the 
Individual  Steel  Bars  (DS  (I) ) 

This  card  is  repeated  for  each  longitudinal  bar 
1-15  Variable  Mean  Value  (F.MEAN  (I )  ) 


1* 

1- 


Format 
F5 . 2 

F  5 . 2 
F 1 0 . 0 
F  1  0 . 0 
F  1 0 . 0 

F5.2 


T4  Variable  Standard  Deviation  (STDV  (I) ) 

31-45  Variable  Constant  (FCONST(I)) 


F5.2 
F  5 . 2 
F5 . 2 
F5.2 
F5.2 
13 

F5 . 2 

F5.2 

t 

F15.5 
F  1 5 . 5 
F 1 5 .  5 


46-50  Variable  Distribution  Type  (ITYPE(I))  15 
Note:  This  card  is  repeated  for  each  variable.  In  this 
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Card 


DATA  INPUT  CONTINUED 

Columns  Data  Description  Format 

program  the  order  of  variables  is  as  follows: 

Concrete  Strength 
Steel  Strength 
Cross  Section  Width 
Cross  Section  Depth 
Core  Width 
Core  Depth 

Distance  From  Compression  Face  to  Nearest 
Longitudinal  Steel 

Distance  From  Compression  Face  to  the 
Longitudinal  Steel  Furthest  From  the 
Compression  Face 

Distance  From  the  Compression  Face  to 
Each  Longitudinal  Bar 


. 

APPENDIX  F 


A  " 

AS 


b" 


^sp 

d” 

D 

e /h 


e 


n 


“ct 


f  • 

c 


'cr 


h 

L 


NOMENCLATURE 

Cross  sectional  area  of  tie  steel,  one  side  of 
column 

Width  of  column  core 

Actual  cover  of  exterior  steel  layers 
Specified  cover  of  exterior  steel  layers 
Depth  of  column  core 
Dead  load 

Eccentricity  of  axial  load  divided  by  the  column 

dimension  perpendicular  to  the  neutral  axis 

Error  in  placement  of  interior  steel  layers 

Modulus  of  elasticity  of  concrete  in  compression 

Modulus  of  elasticity  of  concrete  in  tension 

Modulus  of  elasticity  of  steel 

Concrete  stress 

Concrete  design  strength 

Mean  in-situ  concrete  strength 

Average  concrete  cylinder  strength 

Depth  of  cross  section 

Live  load 


M  Mean  value  of  (B-S) 

(F“  S) 

PACI  ACI  calculated  axial  load 

Ptheory  Axial  load  calculated  from  subroutine  theory 
Ptest  Axial  load  from  Hognestad’s  tests 

K  Nominal  resistance  or  strength 


Spacing  of  ties 


2  07 


208 


NOMENCLATURE  CONTINUED 


£ 

C 


w5oh 


£ 

O 


£ 

t 


£ 

tr 


£ 

u 


5ou 


P 


P 


11 


a 


°(R-S) 


4> 


Coefficient  of  variation  of 
Separation  function,  0.75 
Safety  index 
Dead  load  factor 
Live  load  factor 
Concrete  strain 

Increase  in  strain  at  50%  of  maximum  stress  due  to 

confinement  of  concrete  by  tie  steel 

Concrete  strain  at  maximum  stress 

Concrete  tensile  strain 

Concrete  strain  at  rupture  in  tension 

Crushing  strain  of  unconfined  concrete 

Concrete  strain  at  50%  of  maximum  stress  of 

unconfined  concrete 

Steel  percentage 

Tie  steel  volumetric  ratio 

Standard  deviation 

Standard  deviation  of  (R-S) 

Stress  in  tension 

Rupture  strength  of  concrete 

Understrength  factor 


